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In this article, recent developments in HElium NanoDroplet Isolation (HENDI) spectroscopy are 
reviewed, with an emphasis on the infrared region of the spectrum. We discuss how molecular beam 
spectroscopy and matrix isolation spectroscopy can be usefully combined into a method that provides 
a unique tool to tackle physical and chemical problems which had been outside our experimental 
possibilities. Next, in reviewing the experimental methodology, we present design criteria for droplet 
beam formation and its seeding with the chromophore(s) of interest, followed by a discussion of the 
merits and shortcomings of radiation sources currently used in this type of spectroscopy. In a second, 
more conceptual part of the review, we discuss several HENDI issues which are understood by the 
community to a varied level of depth and precision. In this context, we show first how a superfluid 
helium cluster adopts the symmetry of the molecule or complex seeded in it and discuss the nature 
of the potential well (and its isotropy) that acts on a solute inside a droplet, and of the energy 
levels that arise because of this confinement. Second, we treat the question of the homogeneous 
versus inhomogeneous broadening of the spectral profiles, moving after this to a discussion of the 
rotational dynamics of the molecules and of the surrounding superfluid medium. The change in 
rotational constants from their gas phase values, and their dependence on the angular velocity and 
vibrational quantum number are discussed. Finally, the spectral shifts generated by this very gentle 
matrix are analyzed and shown to be small because of a cancellation between the opposing action of 
the attractive and repulsive parts of the potential of interaction between molecules and their solvent. 
The review concludes with a discussion of three recent applications to (a) the synthesis of far-from- 
equilibrium molecular aggregates that could hardly be prepared in any other way, (b) the study of 
the influence of a simple and rather homogeneous solvent on large amplitude molecular motions, and 
(c) the study of mixed "^He/^He and other highly quantum clusters (e.g., H2 clusters) prepared inside 
helium droplets and interrogated by measuring the IR spectra of molecules embedded in them. In 
spite of the many open questions, we hope to convince the reader that HENDI has a great potential 
for the solution of several problems in modern chemistry and condensed matter physics, and that, 
even more interestingly, this unusual environment has the potential to generate new sets of issues 
which were not in our minds before its introduction. 



I. INTRODUCTION 



A. Scope of the review 



frared spectroscopy has led seamlessly to this apparently 
esoteric, but in reality simple, extension, which, as we 
hope to show here, is likely to spawn a large number of 
scientific and technical applications. 



In this article we review the motivations and the exper- 
imental methodology of HElium NanoDroplet Isolation 
(HENDI) spectroscopy, with an emphasis on the infrared 
region of the spectrum. After reviewing spectrometer de- 
sign criteria and experimental know-how, we summarize 
our present understanding of line spacings, line shapes 
and spectral shifts, and discuss some of the recent re- 
sults in as far as they illustrate general principles and 
possible future applications. Our aim is centered on dis- 
cussing the current understanding of the field, on identi- 
fying unsolved problems and open questions, and not on 
providing a complete summary of the work published so 
far, because the last three years have seen a really explo- 
sive growth in this area. For earlier reviews of this field, 
see Refs. |^, ||. Before describing the motivations and gen- 
eral goals of this kind of research, we shall first show how 
the natural evolution of supersonic molecular beam in- 
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B. The marriage between molecular beams and 
matrix spectroscopy 

Because of the long fluorescence lifetimes of vibra- 
tionally excited molecules, it is impractical, with few ex- 
ceptions ^, to obtain infrared spectra in free jets using 
fluorescence detection. At the same time, when the use of 
well collimated beams is called for, direct measurement 
of the attenuation of the exciting radiation is also out of 
the question 

This leaves bolometric detection of the energy de- 
posited in the beam by a resonant, tunable laser as the 
method of choice for this type of spectroscopy @. As- 
suming unity partition function and full saturation of the 
IR transition, it can be shown that with a seeded beam 
of "standard" intensity and the typical sensitivity of a 
doped Si (or Ge) bolometer operating at 1.5 K (noise 
equivalent power on the order of 100 fW/-\/Hz), signal- 
to-noise (S/N) ratios of up to 10^ can be obtained 

The method is equally applicable to van der Waals 
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complexes which, again with very few exceptions , dis- 
sociate upon vibrational excitation. What is detected in 
this case is a decrease of the total kinetic energy car- 
ried by the beam, because the dissociation fragments are 
deflected from the beam axis, resulting in a decreased 
particle flux to the detector. 

If the complex consists of a large number of monomer 
units, it is not a priori obvious that absorption of energy 
must lead to evaporation, as the amount of excess energy 
per degree of freedom can conceivably be too small. One 
must consider, however, that the temperature of large 
van der Waals clusters is determined by evaporative cool- 
ing ||] which is a self-limiting process. It is then easy to 
appreciate that the absorption of an IR photon by a large 
van der Waals cluster will in general put the cluster above 
the threshold for evaporation to resume. 

Using this detection technique and concepts, IR spec- 
tra of hundreds of van der Waals complexes (mostly be- 
low 4 monomer units) have been obtained during the last 
20 years, contributing a great deal to our knowledge of 
intermolecular forces ^. For larger sizes, many exper- 
iments investigate a single chromophore molecule com- 
plexed to a cluster of rare gas atoms, in order to keep the 
interactions as simple as possible. This class of experi- 
ments has, however, about an order of magnitude lower 
S/N ratios than the experiments carried out on a seeded 
beam of stable molecules. 

In spite of the added difficulties, the possibility of load- 
ing a large rare-gas host cluster with almost any chro- 
mophore by pickup of the latter in a "scattering" cell 
located on the path of the host cluster beam jl^ has gen- 
erated a whole new class of experiments. As the exact 
number of atoms in the cluster starts becoming incon- 
sequential, this approach can be seen as the gas-phase 
analog of Matrix Isolation (MI) spectroscopy |jll[ . As in 
bulk MI spectroscopy, several classes of experiments can 
be carried out that involve both the unperturbed chro- 
mophore and its photodissociation/reaction products. 

In comparison with bulk MI spectroscopy, its Cluster 
Isolation (CI) analogue has two main advantages and two 
major disadvantages. The first advantage is due to the 
finite size of the cluster matrix, which allows control over 
the degree of chromophore complexation. Uncontrolled 
aggregation upon annealing plagues MI spectroscopy and 
sets one of the inevitable limits of this powerful tech- 
nique. The second advantage is the fact that preparing 
the chromophore in a surface location is much easier in a 
cluster than in a matrix. The two disadvantages are (1) 
the usual transient character of every molecular beam, 
which makes it necessary to use powerful laser sources to 
both pump and probe the cluster beam, while the bulk 
matrix can be both manipulated and probed with rela- 
tively low power, broadband light sources, and (2) the 
fact that the temperature of a cluster is "internally set" 
by evaporative cooling and cannot easily be changed. 

Another problem connected with MI spectroscopy is 
the presence of matrix-induced perturbations, which can 
be minimized by the use of bulk solid or liquid helium as 



a medium. This, however, comes with a heavy price: in- 
jection of the sample into helium without causing aggre- 
gation or condensation onto the walls of the container is 
extremely difficult because of low solubility, and has been 
achieved almost exclusively for atoms It should be 
noted, however, that high purity parahydrogen has been 
shown to form stable matrices for molecular dopants with 
perturbations, as measured by widths of rovibrational 
transitions, that are comparable to those observed in he- 
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The extension of CI spectroscopy to helium clusters, 
pioneered at Princeton about a decade ago solves 
the problem of sample injection. It has grown tremen- 
dously in importance since the group of J. P. Toennies in 
Gottingen has shown that a guest in a host helium nano- 
droplet can yield a spectrum indicative of free molecular 
rotation. Almost all the spectra measured to date are 
consistent with the spectral structure predict ed by the 
symmetry of the isolated molecule (see Section [II A ) . In 
contrast, in a solid matrix (or even in a conventional liq- 
uid, if an instantaneous configuration is considered), the 
total symmetry of the system is the combined symme- 
try of the molecule and that of the trapping site, usu- 
ally lower than that of the molecule alone. Further- 
more, the matrix-induced spectral shifts in helium are 
also much smaller than in any other medium, and of- 
ten muc h sm aller than complexation-induced shifts (see 
Section [II E ). For the purpose of comparison with ab 
initio calculations, complexes formed in He droplets can 
then be treated as isolated (gas phase) , and the observed 
shifts provide direct information to decide which isomers 
are actually formed. There are two important implica- 
tions: (1) the structural information contained in the 
spectral symmetry is undiluted, and (2) if a complex is 
formed in helium, the interactions of its parts with the 
surrounding medium are too weak to significantly affect 
its equilibrium structure (zero-point motions and, more 
dramatically, kinetics, are sometimes affected). 

Finally, the possibility of obtaining high resolution 
spectra in helium droplets at 0.38 K and, using •^He, 
0.15 K, allows for the study of superffuidity in nanoscale 
systems, a frontier subject in condensed matter physics 
that has the potential to give a substantial contribution 
to our understanding of highly degenerate many-body 
systems ^ |l|. 



II. EXPERIMENTAL METHODS 

The history of liquid helium cluster/nanodroplet pro- 
duction and the principal methods used for their study 
have been reviewed in the preceding article by Jan 
Northby . Here we will instead discuss the design cri- 
teria of a typical droplet source for HENDI spectroscopy, 
with particular focus on those used in our laboratory. 
The schematic of a typical experimental apparatus is 
shown in Figure |l|. 
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FIG. 1: Schematic of a typical experimental apparatus. From 

Ref. lei. 



A. Nanodroplet beam production 

Due to the weak binding energy of small helium clus- 
ters, large densities and low temperatures are called for in 
the expansion, both resulting in a large gas flux. There- 
fore, liquid helium droplet sources require large pumps 
and/or very small nozzles. The flux is proportional 

to PqcPTq where typical nozzle diameters (d) range 
from 5 to 20 /im, pressures (Pq) from 1 to 100 bar, and 
temperatures (Tq) from 10 to 30 K. Within the ideal 
gas approximation, a 10 /xm nozzle with Tq = 20 K 
and Pq = 15 bar will discharge into the source cham- 
ber approximately the same gas load as a room temper- 
ature 50 /im nozzle at 2 bar, i.e., « 4 standard cm'^/s 
(= 0.16 mmol/s). This is the maximum gas load that can 
be applied to a typical diffusion pump with an effective 
pumping speed of 4000 L/s (referred to N2, i.e., 10000 L/s 
for He) at the limit operating pressure (« 3 x 10~^ torr). 
Typically, droplets with an atom count of a few thou- 
sand are formed. The large mass of these droplets means 
that the beam intensity is not much affected by scattering 
from the background gas. 

If a source is fed at the "standard" load mentioned 
above, it is possible to operate a droplet source with a 
small two-stage closed-cycle refrigerator (the warm stage 
typically operates at 80 K, the cold stage at 12 K) de- 
livering 1 W of cooling power at the cold stage (0.12 W 
of power are needed to cool the gas from 80 K down to 
20 K). However, due to background radiation and back- 
ground gas thermal loads, the use of a 5 W refrigerator 
is recommended. This will allow experiments to span a 
larger temperature and gas load range, and therefore to 
cover a larger range of average cluster sizes. For sources 
which operate below 10 K, which are used for the pro- 
duction of very large clusters ((iV) > 10^), either refrig- 
erators with lower minimum operating temperatures or 
direct liquid He cooling are required. Closed-cycle re- 



frigerators with a zero-load temperature of 3 K are com- 
mercially available of course at substantially higher 
cost. Cost is also the limiting factor for direct liquid He 
cooling, which is the cooling method that allows the max- 
imum flexibility in terms of droplet size range and total 
flux that can be attained ||l9|| . 

In a subcritical expansion, "^He droplets exhibit a log- 
normal size distribution, with a full width at half max- 
imum of 0.55-0.75 (N) 0. While the existence of 
scaling laws for the determination of (N) from expan- 
sion parameters is generally recognized, different formu- 
las have been proposed for it |^ . The most popular one 
states that (N) is a function (not necessarily linear ) 
of Pod^'^T^^^'^. The estimates that we will give here for 
(N) use a formula by Knuth et al. |^^; this formula is 

al- 



in fair agreement with the measurements of Ref. 25 



though the latter seem to scale differently with Tq. The 
case of supercritical expansion, which produces a bimodal 
distribution of much larger droplets, has been analyzed 
in Rcfs. ^ |2^. IR absorption has been used as a 
calibrated source of He evaporation to measure average 
cluster sizes as a function of source conditions p8[ . 

The exact shape (density distribution) of a He droplet 
has been computed using both density functional [^0|, |2^ 
and Monte Carlo ||2^ methods. A good approxima- 
tion for "back of the envelope" calculations is that of 
a sphere having the same density as bulk liquid helium 
(po = 0.0218 k-^, for -^He [|o|). According to this liquid 
drop model, a droplet containing TV '*He atoms has a ra- 
dius R = 2.22iVi/3 A. A more realistic description uses 
the analytic function pG, pi[: 



p[r) 
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which accounts for the diffuseness of the surface via a 
thickness parameter g; the 10-90% thickness (« l.lg j is 
estimated by various methods to be « 6 A for "^Hc EG] 
and « 8 A for ^He 

The equilibrium temperature of a droplet is dictated 
by evaporative cooling and can be extracted 

from the fit of the ro-vibrational spectrum of a probe 
molecule (Sec. II Dl). The measured values, 0.38 K for 
^He (Ref. 34, and most entries in Tab. [|) and 0.15 K for 
■^He [ p4[, are in good agreement with theoretical predic- 
tions ]32| , |3^ . 

About 1 cm downstream from the nozzle, the droplet 
beam is collimated by a skimmer. As shock waves may 
play a role in attenuating the beam, it is important that 
the shape of the skimmer is optimized, and that its edges 
are as sharp as possible. Typically, commercial electro- 
formed conical skimmers, about 500 /xm in diameter, are 
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To give an idea of what intensity can be expected for a 
droplet beam, consider a detector located 500 mm away 
from the source, with a 1 mm^ detection area, and an 
average cluster size of {N) 4000. Assuming that 10% 
of the flux is in the droplets and that the droplet flux is 
spread uniformly over 1/10 of a steradian in the forward 
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direction (a conservative estimate) , one gets a flux of on 
the order of 10^^ droplets per second at the detector 
Notorious difficulties in calibrating the absolute values of 
molecular beam detectors (be it bolometers or mass spec- 
trometers) make this number hard to calculate precisely. 
For larger clusters, a first approximation of the observ- 
able mass flux can be obtained by scaling the number of 
clusters inversely with the number of atoms per cluster. 

Whereas pure and doped "^He droplets have been ex- 
tensively studied far fewer experiments with pure '^He 
or mixed ^He/^Ile droplets have been done |2|, 
H pL E^' El chiefly due to the high cost of 
the isotope ^2[. Budget dictates that •^He be continu- 
ously recycled, requiring specialized equipment: sealed 
pumps and a purification system. The gas is collected 
from the exhaust of the source chamber pumping sys- 
tem, cleaned on liquid-nitrogen-cooled zeolite traps, pres- 
surized by a membrane compressor to 20-80 bar and 
fed back to the nozzle. The recycling system used in 
Refs. m, H, H, 1^, im needed about 20 bar-L gas to 
be filled. The loss rate under normal operation was < 2 
bar-L per week Thus although a considerable initial 
investment is required the operating costs are comparable 
to the other running costs of these experiments. 

There is an important difference between '^Hgn and 
■^Hbat clusters: whereas the former are bound for any 
A^lE^, calculations predict that small •^He droplets of 
less than about 30 atoms are unstable Q, ^ because of 
the larger zero point energy associated with the lighter 
isotopic mass. Experimentally, a strikingly different be- 
havior of average droplet size as a function of nozzle 
temperature is observed for the two isotopes. Whereas 
the size of "^He droplets increases gradually with de- 
creasing temperature the formation of "^He droplets 
sets in sharply, and at a considerably lower temperature 
(To = 11.5 K for Pq — 20 bar and d — 5/xm, giving 
droplets of about 4500 atoms ^^)- In the expansion 
of an isotopic mixture, formation of *He-enriched mixed 
droplets is observed The ^He concentration in 

the droplets is several times higher than that of the ex- 
panding mixture (< 4%) and depends on the expansion 
conditions . Because ^He will selectively solvate most 
dopants, the production of pure '^He droplets requires gas 
of high isotopic purity ("^He < 10~^), which is commer- 
cially available M]. 



B. Nanodroplet beam doping 

After coUimation, the droplets are seeded with the 
molecule of interest by passing them through a small gas 
pickup cell where the pressure is tailored to obtain the 
desired mean number of dopants per droplet. The in- 
ternal and kinetic energy of the dopant, as well as the 
energy of solvation, are dissipated into the droplet and 
result in evaporation of He atoms (200-250 atoms, each 
taking off 5.5-7 K of energy |Q, are estimated for the 
pickup of a small molecule such as HCN), leading to a 



microcanonical system cooled to 0.38 K (we assume pure 
■^He droplets here). On average, the impact parameter of 
the molecule is significant (2/3 of the droplet radius), re- 
sulting in a mean angular momentum of « 45iV^/'^?i per 
collision with HCN. Likewise, He evaporating at 0.38 K 
will carry away 10-15?i of angular momentum per atom. 
At present, it is not clear how much angular momentum 
remains in the cluster after evaporative cooling, nor in 
which form it is stored if it exceeds the thermal values. 

A quick estimate of the optimum pickup conditions 
can be made by assuming unity sticking coefficien t an d 
a droplet cross section of 15.5A^^/"^ (see Section |H A| ); 
further, the droplet's size reduction upon pickup and the 
nonzero velocity of the gas molecules are neglected (if 
the latter is account ed for, the co llision rate will increase 
by a factor of « y/l + (i^m/wd)^, where Vd « 400 m/s 
is the droplet velocity and w,„ is the rms velocity of the 
molecules in the pickup cell) . With these approximations 
(further corrections are discussed in Ref. 47), the pickup 
process results in a Poisson distribution of the number 
of dopant molecules per droplet, and a column density 
corresponding to « 2.57V~^^'^ Pa-cm is required to maxi- 
mize the probability of single pickup. Under typical con- 
ditions, the vapor pressure needed for efficient doping of 
helium droplets is on the order of 10^^ Pa, some four 
orders of magnitude lower than what is typically used 
in supersonic jet coexpansions. This implies that many 
species can be doped into helium that are not practical 
to study by jet spectroscopy because of the high tem- 
peratures involved, which are either difficult to reach or 
cause thermal decomposition of the species. The large 
cross sections for the droplets imply that the vacuum re- 
quirements in the instruments are more stringent than 
for beam experiments on monomers or small clusters. 

Molecules captured by the same droplet are expected 
to find each other and form van der Waals complexes, but 
once past the pickup region, each droplet becomes an iso- 
lated system and no further aggregation occurs. At the 
temperature of the droplet, 0.38 K, the ro-vibrational 
spectrum of a typical molecule spans ~ 1 cm~^, which 
is generally smaller than the shifts induced by complex- 
ation. Hence, spectra of different small oligomers will 
not overlap, and pickup of more than one molecule per 
droplet will only decrease the intensity of the single- 
molecule spectrum, without affecting its shape. For this 
reason, at least when one is concerned with monomers 
and dimers only, high purity of the dopant gas is de- 
sirable but not essential. An extreme case is the spec- 
trum of H^^C-^^CH, which has been observed Q at 
its naturally occurring abundance (2.2%) interleaved to 
that of H^^C^^CH. The situation is different for higher 
oligomers, because complexation shifts tend towards an 
asymptotic limit |5^, and the spectra tend to merge. 
Nevertheless, in a favorable case complexes of OCS 
with up to 17 hydrogen molecules have been made and 
spectroscopically resolved though this was done in 
mixed '^He/^'He droplets whi ch coo l to 0.15 K. Pendu- 
lar spectroscopy (see Section IID4 and Ref. Eol) allows 
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even the modest width of the rotational contours to be 
collapsed, which is particularly favorable in separating 
closely spaced oligomer or isomer bands. 



C. Radiation sources 

Since spectroscopy of doped helium droplets is still a 
relatively new field, it is often the availability of a suitable 
radiation source that dictates the choice of the system to 
be studied. In the early days SFg {1^3 mode) 

was chosen as a dopant molecule because of its large ab- 
sorption cross section in a spectral region where watts of 
radiation (from line-tunable CO2 lasers) were available. 

To achieve rotational resolution, this transition was 
later studied with a continuously tunable semiconductor 
(lead salt) laser |Q. The same type of laser was also em- 
ployed for the beautiful studies of OCS in mixed ■^Hc/^'He 
carried out in Gottingen |15, ^ |4[ p5| . 

Pulsed lasers have also been used in early experiments 
on HF-, H2O- or NHa-doped helium droplets fs^, |8). 
They offer the advantage of high power and wide-range 
tunability; those come at the price of usually broadband 
emission, so that it becomes difficult to establish whether 
the width of the observed spectral lines is determined 
by the surrounding helium environment or rather by the 
laser itself (via instrumental linewidth or power broad- 
ening) . 

High power line-tunable gas lasers, cw or pulsed, are 
still being used up to the present time ^ , 61 , 32j . Be- 
cause of the sparse set of available lasing frequencies, they 
are most suited to the investigation of small molecules, 
with small moments of inertia and therefore well spaced 
lines. Thanks to their large output power, they have 
been successfully used to perform saturation measure- 
ments ||2|. 

Most of the experimental spectra to date have been 
acquired with color center lasers, in the 3 /xm re- 
gion [|9[ 1^, |6|, H, MM, @ ' and, to a lesser extent, 
in the 1.5 /im region]|69|T^ These lasers offer broader tun- 
ability range, comparatively large output power (30 and 
250 mW, respectively) and high resolution, as they can 
be stabilized to better than 1 MHz without too much 
difficulty. Since the 1.5 /Ltm region is used in telecommu- 
nications, commercial optical components (e.g., fibers, 
connectors, and high finesse cavities) are available for 
this type of measurements. 

At these power levels it is possible to bring near sat- 
uration a 10 MHz homogeneously broad, strong funda- 
mental transition near 3 /xm. If the lines are broader 
or the transition weaker (as in the case of overtones at 
1.5 fiTo), laser field enhancement devices can be usefully 
employed (see Section II D). Looking at the future, since 
HENDI spec tral li nes are rarely narrower than 0.01 cm~^ 
(see Section [II C ) and since pulsed lasers with similar 
bandwidth are available, it would probably be useful to 
explore the possibility of developing pulsed helium nano- 
droplet sources. Pulsed devices that produce small He 



clusters are already available at present |70[ . 

For cw beam sources, narrow bandwidth high repe- 
tition rate pulsed lasers have recently become available 
at power levels that allow frequency multiplication and 
mixing. As the use of regenerative amplifiers has allowed 
repetition rates up to 10 kHz, the problem of the low duty 
cycle characteristic of traditional pulsed laser sources can 
be overcome. For strictly cw laser sources, recent devel- 
opments that will likely find HENDI applications are cw 
OPOs ||7l| and quantum cascade lasers ||7^. Quantum 
cascade lasers offer the advantage of being available to 
wavelengths as high as « 24 /xm , which is going to be 
particularly useful for the study of large molecules where, 
at higher frequencies. Intramolecular Vibrational Redis- 
tribution (IVR) is an unavoidable source of line broad- 
ening and of loss of information. 



D. Excitation schemes 



1. Infrared 
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FIG. 2: Spectrum of the 21/1 transition of propyne, and fit to 
a symmetric top. The inset shows the hnewidth for each tran- 
sition. Note that at the temperature of the droplet (0.38 K) 
the Q-branch would be absent if nuclear spin relaxation had 
taken place. From Ref. 

Infrared (IR) spectroscopy in He droplets derives its 
uniqueness from the fact that the medium still allows ro- 
tational resolution (a representative spectrum is shown 
in Figure |^) , even for relatively large molecules and com- 
plexes. This is to be compared to "classical" solvents, 
for which rotational resolution is, at best, only possible 
for molecules of the smallest moments of inertia. Rota- 
tional resolution allows the measurement of the cluster 
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temperature, a very sought after quantity as the follow- 
ing statement (made in 1994) proves: 

But so far no "thermometer" to measure a 
cluster's temperature is available. Worse, 
there does not even exist a suggestion how 
to construct one. 

[Ref. Chapter 1.5: Experiments not possi- 
ble today] 

Ironically, in the same book, the first results on infrared 
spectroscopy in He clusters were summarized, and the 
following prediction was made: 

In future the technique may be combined 
with continuously tunable narrow band lasers 
of higher power as these become available. 
[ibid., Chapter 3.7: Infrared spectroscopy] 

Continuously tunable, narrow-band lasers were indeed 
the key that, one year after, allowed for rotational reso- 
lution to be achieved and thus temperature to be mea- 
sured |5^, The technical solution that was used 
to overcome the power requirement consists of a differ- 
ent arrangement of how the laser beam interacts with 
the beam of doped droplets. In contrast to the opaque 
bolometer detectors of the early experiments, the use of 
a "transparent" mass spectrometer detector allowed for 
a setup where the laser and cluster beam counterpropa- 
gated, resulting in interaction region of tens of centime- 
ters. With this arrangement, low power diode lasers can 
be profitably employed . Additional benefits of using 
this setup for pulsed lasers is that their duty cycle (de- 
fined as the fraction of the beam that is in the laser field) 
is increased from 10~^ to 10"'^ and the detector output 
can be time gated to match the depletion transient. 

When bolometers are used as detectors, multipass 
cells Q and power build-up cavities have been suc- 
cessfully employed. While a well aligned multipass cell 
can give a 30- to 50-fold increase in signal relative to 
the single-crossing case, a resonant cavity can further en- 
hance the radiation intensity seen by the droplets by one 
order of magnitude. Naturally, this comes at a price, 
since, as the laser is scanned, the cavity has to be kept in 
resonance by suitable feedback electronics. Moreover, as 
currently implemented, the interior of the cavity is a vol- 
ume bound by metallic surfaces; as such it is effectively 
shielded from electromagnetic fields. Hence, double res- 
onance and Stark spectroscopy experiments are not pos- 
sible. Nevertheless, the large amount of available power 
(currently only matched by line-tunable gas lasers j62|) 
makes this approach very desirable to reach saturation 
conditions, and to study forbidden transitions |69| . 

HENDI spectra have been measured by beam deple- 
tion using both bolometers and mass spectrometers to 
detect the on axis beam flux. Even with helium evapora- 
tion occurring at 0.38 K, the mean lab frame scattering 
angle of the evaporated helium is still « 10 degrees, and 
so the evaporated helium largely miss the detector. The 
relative sensitivity of the two detection methods is hard 



to determine unambiguously from the literature spectra, 
as they involve widely different source powers, transition 
moments, and detection geometries. Our attempt to cor- 
rect for these factors suggests that detection of depletion 
with a bolometer is about a factor of ten more sensitive 
than with a mass spectrometer. This assumes use of a 
continuous wave laser and of a multipass cell. The lat- 
ter partially offsets the shorter excitation length (relative 
to the more efficient counter propagation scheme possi- 
ble with the mass spectrometer) so that ultimately only 
modest reduction results. The mass spectrometer will 
likely give superior sensitivity in the case when one is us- 
ing a low repetition rate pulsed laser for excitation, due 
to the possibility of using gated detection. 



2. Microwaves 

Investigation of purely rotational transitions in doped 
helium clusters shifts the frequency spectrum to the mi- 
crowave (MW) region. Radiation is confined inside a 
section of waveguide (part of which is coUinear with the 
cluster beam) properly terminated at the ends. Confine- 
ment is necessary not only to concentrate the available 
microwave power, but also to minimize its absorption or 
rectification by the bolometric detector, which can easily 
obscure the true signal. 

Pure rotational spectroscopy in He droplets differs 
from vibrational and electronic spectroscopies in that it 
is inherently multiphoton. One can see from Table ^ that 
in a He droplet the energies associated with a pure ro- 
tational transition (0.5-70 GHz) are insufficient to evap- 
orate even a single He atom. Microwave spectroscopy 
in a He droplet relies on the occurrence of rotational re- 
laxation in an isolated droplet, which allows a dopant 
molecule to repeatedly absorb photons and release their 
energy into the droplet. The fact that the observed sig- 
nals are of the same magnitude as for infrared transitions 
has been used to estimate the number of relaxation cy- 
cles that each molecule must undergo in the « 100 /is 
it spends in the MW field, and thu s to place the ro- 
tational relaxation time (see Section [ HI Cl ) on the sub- 
microsecond time scale for HCCCN |76| . 

The large amount of available microwave power (tens 
of watts), along with the large transition dipoles in this 
spectral region has made saturation measurements pos- 
sible ^ , which have been used to establish and 
quantify the dominance of inhomogeneous broadening in 
HENDI spectroscopy. 



3. Double resonance 

The term double resonance (DR) refers to the cou- 
pling of two transitions by a common level, which al- 
lows driving one transition (the "pump") to change the 
strength of the other (the "probe"). In the gas phase, 
DR is routinely used (1) to simplify the assignment of 
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complex spectra, (2) to reach states inaccessible by one- 
photon selection rules, and (3) to measure homogeneous 
lines within an inhomogeneous profile. To date, with one 
exception, only type (3) experiments have been done in 
helium droplets, and only in the simplest "steady state" 
versions, i.e., without explicit time resolution of the time 
delay between pump and probe fields. The exception is 
the recent hole burning study of the electronic spectrum 
of tetracene 78 which shows a helium induced splitting 
of the electronic band origin. The complexity of the dy- 
namics in helium droplets has unfortunately made the 
interpretation of these experiments much more difficult 
than that of their gas phase analogues. 

The first double resonance study in helium were the 
MW-MW experiments performed on the R(3) and R(4) 
transitions of HCCCN |Q . Saturation measurements on 
these rotational lines had shown the typical pattern of a 
line to be dominated by inhomogeneous broadening. It 
was expected that DR spectra would directly reveal the 
homogeneous component of the line, as in a classical hole- 
burning experiment (Ref. page 438). The observation 
of lines much broader than expected was interpreted as 
a sign of spectral diffusion (see Section [II C 2). 

In an attempt to gain further insight, MW-IR exper- 
iments on HCCCN [Q have been performed. Indepen- 
dent of this work, MW-IR double resonance experiments 
where performed on OCS Because rotational relax- 
ation times in He clusters lie in the nanosecond range 
(during which a cluster only travels a few /im), double 
resonance experiments require spatial overlap of the two 
radiation fields. While in MW-MW experiments this is 
readily achieved, in MW-IR experiments spatial overlap 
has been achieved by either propagating the laser along 
the axis of the waveguide Q, or by using the walls of 
the waveguide as a substrate for thin mirrors in a mul- 
tipass arrangement (see Figure p. The long wave- 
length of microwaves, compared to the transverse dimen- 
sions of the cluster beam, makes it possible to have holes 
drilled into the waveguide as ports for the cluster and 
laser beam, without introducing significant microwave 
losses. 

While a MW-IR scheme introduces the additional 
complication of vibrational excitation (which fans out 
the "probe" transition frequencies associated to a given 
"pump" frequency), it has the advantage that the probe 
laser can cover the entire manifold of rotational levels, 
thus probing rotational states that are not directly con- 
nected via the pump field. With some differences, prob- 
ably due to the use of different molecules, both exper- 
iments find that the population of all levels is to some 
extent affected by the pump beam, which indicates fast 
rotational diffusion. 



4- Stark spectroscopy 

Stark spectroscopy has been shown to be a potentially 
valuable tool to learn more about the interactions of the 



dopant molecule with the helium environment and possi- 
ble distortions of the molecule (particularly of complexes) 
by gentle helium solvation ^ . In a simple linear di- 
electric approach, one would expect the molecular dipole 
moments to be reduced in helium by its dielectric con- 
stant, 1.055 However, this is only valid for a dipole 
in a continuous medium. If the dipole is located in the 
center of a spherical cavity, the induced dipoles in the 
solvent cancel out and give no net change in dipole mo- 
ment. As the solvation around molecules is in general 
highly anisotropic, the net induced moment provides a 
potentially interesting measure of the anisotropy. Re- 
cently such calculations have been performed for Mg (3P) 
solvated in He, with the anisotropic density obtained by 
a Density Functional method, and a 6.2% reduction of 
the optical transition dipole of Mg (3P — > 3S) transition 
has been calculated [g2|. Such a reduction, along with 
the solvent-induced red shift, has been found to quantita- 
tively account for the experimentally measured increased 
lifetime of the transition. The solvent-induced stabiliza- 
tion of dipolar structures was listed as one of the possible 
causes of the decreased tunneling splittin g obs erved for 
the HF dimer in He droplets (see Section IV E| ). 

Nauta and Miller have studied the Stark pattern 
of the R(0) line of HCN in the IR. They find the split- 
ting and intensity pattern expected from the isolated 
molecule, but with the observation that what appears to 
be the a splitting of the \M\ = 0, 1 components persists 
even at zero field, as a poorly resolved doublet structure 
of the R(0) fine (see Figure ||). A careful study of the 
dipole moment in this and other systems would likely be 
quite valuable for testing theoretical predictions of the 
helium solvation structure. 

Another relevant manifestation of Stark effects is the 
mixing of rotational states (Ref. ^ Section 10.6). In 
He droplets, it has been exploited to induce a (normally 
forbidden) Q-branch in those linear molecules, such as 
HCN and HF, whose zero-field rotational constants are 
large enough that their spectrum consists of the R(0) line 
only. By this method, the band origin and the rotational 
constant of the molecule can be directly measured Q . 
In this application the field is kept as weak as possible, 
compatibly with the need of inducing a measurable sig- 
nal, and indeed an accurate value for the band center is 
obtained by measuring the position of the Q-branch at 
increasingly weak fields and extrapolating to zero field. 

The electric field can be used to "tune" the frequency 
of a given transition. This has been exploited to quench 
tunneling in (HF)2 p5[ , and to perform time-resolved 
relaxation measurements in the ms range (by mapping 
different excitation times onto an electric field gradi- 
ent) [H. 

Once the information contained in the ro-vibrational 
Hamiltonian has been extracted, rotational resolution ac- 
tually becomes a nuisance: first, the band strength is 
fractionated, resulting in smaller signals; second, a typi- 
cal spectrum spans a range comparable to fceT (0.38 K 
or 0.26 cm""'^), which results in a corresponding loss of 
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resolution when trying to identify closely spaced bands (a 
common occurrence when many different van der Waals 
complexes of the same molecules are present). For gas 
phase polar molecules, an elegant solution is Pcndular 
State spectroscopy fss] , |8^ , which is the application of a 
static electric field sufficiently strong to convert the rota- 
tion of a molecule into oscillations of the dipole moment 
around the external field. In the strong static field limit 
(and with the permanent and transition moments par- 
allel to each other), the entire rotational contour of the 
spectrum collapses into a single line at the wavenumber 
of the purely vibrational excitation, with an integrated 
intensity 3 times that of the entire band at zero field 
(the traditional spectroscopic "conservation of intensity" 
holds, it is just the orientation of the molecules is now 
fixed, and absorption, no longer averaged, will be maxi- 
mum when the transition dipole and the laser polariza- 
tion are parallel). It is indeed possible to use the signal 
strength as a function of relative polarization of the op- 
tical and Stark fields to determine the angle between the 
static and transition dipoles. 

This technique has been ported to He droplets by 
Nauta and Miller It is particularly attractive 

in this environment because both the low temperature 
and the increased moments of inertia make it easier to 
reach the required "high field" limit. By this technique, 
the spectra of linear chains of up to 8 HCN units p9[ , 
and up to 12 HCCCN units Q have been resolved (see 
Section IV A and Figure |). 
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FIG. 3: Stark spectrum of linear chains of cyanoacetylene 
(labeled by the number of monomer units), showing the high 
resolution achievable with pendular state spectroscopy. The 
series is taken at progressively smaller droplet sizes, showing 
that the size of the droplet limits the length of the chains that 
can be formed. From Ref. 

Suppression of rotation can also be achieved by re- 
placing '^He with '^He |5^. The technique has been 
used to collapse the absorption bands of 0CS-(H2)„ com- 
plexes so that n could be assigned via the measured band 
shift ^ . This method is applicable even to nonpolar 
cases, but lines will typically not be as narrow as for "^He 



spectra. 



III. CURRENT UNDERSTANDING AND OPEN 
QUESTIONS 

A. Molecular symmetries in the nanodroplet 
environment 

The most remarkable feature of the vibration-rotation 
spectra of molecules in liquid ^He nanodroplets is the 
presence of rotational fine structure, either quantum 
state resolved or observed as rotational contours. This 
was first seen in the seminal work of Hartmann et al. 
on the spectrum of SFg, but has now been observed a 
wide range of other molecular solutes. Table | contains 
a list of molecules studied to date with rotational resolu- 
tion in the IR, at least as known to us. In almost every 
case, the nature of the spectroscopic ro- vibrational struc- 
ture is exactly what one expects from the gas phase, i.e., 
the standard spectroscopic signatures of a linear, sym- 
metric, or spherical top are preserved, along with the ro- 
tational selection rules based upon the symmetry of the 
gas phase transition dipole moment, though often with 
dra matica lly changed spectroscopic constants (see Sec- 
and Figure |^). High temperature spin statis- 
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tical weights have been observed, indicative of negligible 
relaxation between nuclear spin isomers on the time scale 
of the experiment (~ 100/xs), as first noted by Harms et 
al. Q (see also Figure ||). 

Viewed in terms of a classical picture of the instan- 
taneous configurations of the helium atoms, the above 
preservation of symmetry is unexpected in that these 
instantaneous configurations do not retain the molec- 
ular symmetry. In view of the large delocalization of 
the He atoms, we must rather consider the symmetry 
of their ground state many-body wavefunction, which is 
the same as the molecular symmetry when the helium 
is in its ground state. Bulk excitations of the helium 
droplet would lower this symmetry, but at the temper- 
ature of the droplets, liquid-drop-model calculations in- 
dicate that there are no thermally excited phonons [ p3[ . 
Thermal excitations are present in the surface modes of 
the droplet, but these appear to couple weakly to the 
solvated molecules, though they must ultimately be re- 
sponsible for the relaxation of the latter into thermal 
equilibrium with the droplet. Thus, we are lead to view 
the helium not as a classical fluid of rapidly moving par- 
ticles, but as a quantum probability density associated 
with a bosonic iV-body wavefunction, which is a smooth 
and continuous function of coordinates and reflects the 
symmetry of the potential that the molecule creates for 
the helium. 

One potential counterargument of the above model is 
that linear molecules will create a cylindrically symmet- 
ric ring of helium density around them, implying that 
the entire system should act as a symmetric top, which 
has an additional rotational degree of freedom. However, 
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it is easily shown that the Bose symmetry of the hehum 



TABLE I: Rotational constants derived from rotationally re- 
solved IR or MW spectra. Blank entries in the "Mode" col- 
umn indicate purely rotational spectra. The modes in paren- 
thesis indicate that more than one band has been seen but the 
measured rotational constants do not differ enough to grant 
separate entries. 



Molecule 



DCN 
HCN 

HCN 
(HCN)2 
(HCN)3 
HCCCN 
HCCCN 
(HCCCN)2 (A) 
(HCCCN)2 (B) 
HCCCCH 

jj13q13qjj 

CH3CCH 
CH3CCH 
CF3CCH 
CF3CCH 
HCN-N2 

HF 

(HF)2 
CO-HF 
N2-HF 
OC^^S 

oc^^s 

pHa-OC^^S (A) 
pH2-OC^''S (B) 
pH2-0C^''S (C) 
SFe 

HCOOH (A) 

HCOOH (B) 

HCOOH (C) 

DCOOH (A) 

DCOOH (B) 

DCOOH (C) 

CO2 

NHa 

0H2-HF 

pDz-HF 

0D2-HF 

CH4 

N2O 

(CH3)3SiCCH (A) 
(CH3)3SiCCH (B) 
Ar-HF 

Cyclopropane {A + 
Cyclopropane (C) 
0H2-HCN 
pDz-HCN 
Mg-HCN 
Mg3-HCN (A) 
Mg3-HCN (B) 



Mode 


^gas 




cm~^ 




1.208 




1.478 




1.478 




0.0582 




0.0156 




0.1517 












0.0113 


V3 


0.1464 


V3 


1.1766 


V3 


1.1195 


Vl 


0.2851 


2vi 
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Buc 
cm~ 



Ref. 



2z^i (A) 
2vx (B) 

1^1 {1^2) 

Vl 
Vl 
V3 
V3 
V3 
V3 
V3 
V3 
V2 {V\) 



0.1908 

0.0960 

0.0525 

19.787 

0.2167 

0.524 

0.1066 

0.2029 

0.1979 



0.0911 
2.5758 
0.4020 
0.3470 



V3 

V2 
V\ 
V\ 
V\ 
V3 

2ui_ 



0.39 

9.96 

0.8325 

0.492 

0.487 

5.2406 

0.419 

0.1056 

0.0654 

0.102 

0.6702 

0.4188 

0.4303 

0.2663 



0.999 

1.204 

1.175 

0.0193 

0.0073 

0.0525 

0.0518 

0.0168 

0.0050 

0.047 

1.0422 

0.987 

0.0741 

0.0717 

0.1004 

0.0355 

0.019 

19.47 

0.0986 

0.1022 

0.045 

0.0732 

0.0706 

0.0847 

0.0544 

0.0422 

0.034 

1.3811 

0.2970 

0.1996 

1.16 

0.24 

0.21 

0.154 

7.5 

0.182 

0.164 

0.159 

5.013 

0.0717 

0.0275 

0.0144 

0.04 

0.1631 

0.2416 

0.0897 

0.0785 

0.0285 

0.035 

0.0167 
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eliminates the low lying K excited states : states with 
angular momentum along the symmetry axis cannot be 
thermally populated, and the K = Q spectrum of a sym- 
metric top is indistinguishable from that of a S state 
of a linear molecule (see also Section IV D). This argu- 
ment alone does not exclude the possibility of creating 
K > Q states with modest energy (~ 3.2 K has been es- 
timated for OCS in He |]88| ). However, in order to create 
such a "vortex" excitation localized in the first solvation 
shell, it is necessary to have a nodal surface separating 
i t fro m the vortex-free remainder of the helium droplet 
[149 1; such a nodal surface may come at a steep energy 
cost, making the creation of high-if modes even more 
difficult. It would be interesting to have an estimate for 
the excitation energy of such states calculated with the 
fixed node diffusion Monte-Carlo (DMC) method. 

Roger Miller's group has observed two examples of 
spectra of highly polar linear molecules (cyanoacetylene 
dimer and HCN trimer [50[|) which appear to have Q- 
branches, and thus suggest a breakdown of this "rule" 
of preservation of molecular symmetry. The origin of the 
Q-branches in these cases has not yet been explained con- 
clusively. It has been proposed that they arise not from 
population of higher K states of the helium, but are due 
to some mechanism that would create an anisotropic po- 
tential for rotation of the molecules. One candidate is the 
anisotropy due to the finite size of the droplet and the 
delo cahza tion of the molecule in the droplet Q (see Sec- 
This model predicted the wrong droplet size 
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dependence of the Q-branch intensity unless an ad hoc 
energy term was added to favor location of the molecule 
slightly below the helium surface. Another proposal 
is the presence of a vortex in a minority of droplets, which 
is expected to strongly align molecules along it |Q . This 
would also appear to favor Q-branch intensity in smaller 
droplets since a larger fraction of pickup collisions in 
such droplets would be expected to impart enough an- 
gular momentum to create a vortex (on the order of one 
fi per helium atom). We would like to add to the list 
the anisotropy that arises from the hydrodynamic cou- 
pling of rotation and translation |9^ (see Section III E| ), 
which is expected to be most important for precisely the 
extremely prolate molecules for which this interesting ef- 
fect has been observed. Furthermore, this effect should 
be most important for large droplets, where the motional 
averaging of the linear momentum (to which the angular 
momentum is coupled) is slowest. Clearly, these "odd" 
spectra should be subjected to further experimental and 
theoretical work since some important new dynamical ef- 
fect may be revealing itself by these phenomena. 

It is worth noting that none of the ro- vibrational spec- 
tra observed in helium have shown evidence of "phonon 
wings" . Such features, which appear to the blue of the 
molecular excitation and involve creation of one or more 
helium phonons upon solute excitation, are observed 
and often predominate in electronic spectra p3[ . 
The strength of the pure molecular excitation, the zero 
phonon line (ZPL), is expected to decrease exponentially 
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with the degree of solvent reorganization upon molecu- 
lar excitation, and thus the complete dominance of the 
ZPL in the case of vibrational excitation implies that the 
reorganization of the solvation density must be negligi- 
ble relative to the zero point fluctuations of the helium 
ground state density. Perhaps the study of flo ppy, large 
amplitude vibrational modes (see Section tVB) will show 
evidence of the phonon wings, which could provide valu- 
able data on solvent reorganization in a quantum liquid. 



B. Translational motion of molecules: Where are 
the molecules? 

The gas-phase-like symmetry of ro-vibrational spectra 
of molecules attached to helium droplets strongly sug- 
gests that these molecules are located in a nearly isotropic 
environment. This rules out surface sites since their re- 
duced symmetry would lead to a characteristic splitting 
of the M degeneracy, as has been seen for the electronic 
spectrum of He2 which resides above the surface of the 
helium jo^. Once a molecule penetrates the helium sur- 
face by more than the first or perhaps second solvation 
layer, a relatively flat potential is expected. Vilesov and 
Toennies proposed treating the motion of the sol- 
vated molecule as that of a free particle in a spherical 
box. However, one can explicitly calculate the effective 
potential that pulls the molecule toward the center of the 
droplet by summing up the missing long range interac- 
tions outside of the droplet (i.e., taking the molecule in 
bulk helium as the reference state). The resulting po- 
tential turns out to localize molecules well within the 
droplets (see Figure ^ , keeping them away from the sur- 
face region where the approximations used to derive the 
potential are expected to break down ]9l| | . The spectrum 
of thermally populated energy levels of a particle in this 
potential resembles more closely that of a 3-dimensional 
harmonic oscillator than that of a particle in a box; the 
heat capacity arising from this motion is nearly Sfce and 
is similar to that of an entire, pure helium nanodroplet of 
3 nm radius. The effective frequency for vibrational mo- 
tion in this potential is « 0.5-1 GHz in typical cases. The 
potential is proportional to the molecule-helium Cg coef- 
ficient minus the helium-helium Cg coefficient times the 
net number of helium atoms displaced by the molecule, 
i.e., there is a buoyancy correction due to the change 
in the net helium-helium interaction as the molecule is 
displaced p6|| . 

For a non-spherically-symmetric molecule, displace- 
ment from the center of the droplet leads to a reduction in 
site symmetry and a weak coupling between translational 
and rotational motion. The leading coupling term has a 
P2 (cos 9) dependence on the angle 9 between the axis of 
a linear rotor molecule and the displacement vector from 
the center of the droplet (P„ are the Legendre polyno- 
mials). If the displacements from the center are assumed 
to be static and their distribution given by Boltzmann 
weights, this P2 asymmetric field results in the same type 



of characteristic splitting of M components of a tran- 
sition as discussed above for a surface location, and is 
easily calculated analytically. However, calculations in- 
dicate that the expected splittings are on the order of, 
or less than, the vibrational frequency in the 3-D poten- 
tial, and thus the splittings are predicted to be at least 
partially motionally averaged. Lehmann | pl| has given 
explicit matrix elements for calculations in a coupled ba- 
sis of translational and rotational angular momenta of 
the molecule, and presented some sample calculations. 

Another important coupling of translational motion 
and molecular rotation is exp ected to arise from a hydro- 
dynamic effects (see Section HI D 1 ) . The translational 
motion of an ellipsoid in helium leads to a hydrodynamic 
contribution to the effective mass, which depends upon 
the orientation of the ellipsoid relative to the linear mo- 
mentum. Lehmann | ]9l[ | has derived the quantum Hamil- 
tonian including this coupling term, and the resulting 
matrix elements between the translational and rotational 
angular momenta. Model calculations on simple systems, 
such as HCN and OCS, suggest that the hydrodynamic 
coupling is more important than the potential terms men- 
tioned in the previous paragraph. It is interesting to note 
that the hydrodynamic term leads again to the expected 
P2 (cos 9) induced splittings in the limit of bulk helium, 
where the potential asymmetry goes to zero. 

The relatively dense manifold of center-of-mass states, 
coupled to the rotation of the molecule which changes 
upon ro-vibrational excitation, creates a source of inho- 
mogeneous broadening in the spectrum that was origi- 
nally not recognized: quantization of the translational 
motion should lead to a potentially resolvable fine struc- 
ture of the spectral features. However, it is likely that 
the broad droplet size distribution present in all ro- 
vibrational spectra observed to date will wash out any 
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FIG. 4: Isotropic (V2) and anisotropic (V4) parts of the con- 
fining potential for a dopant molecule (HCN) in He droplets, 
arising from the dopant-He long-range interaction. P{r) is 
the radial probability density function for the position of the 
dopant. From Ref. bll. 
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details of this fine structure, which is then only observ- 
able as a source of inhomogeneous broadening in the spec- 
trum. 

These couplings of the center of mass and rotational 
motion of the molecules are believed to make important 
contributions to the inhomogeneous lincshapcs often ob- 
served in rovibrational spectra. However, calculations 
done to date have neglected the vibrational dependence 
of the interaction of the helium with the molecules. Any 
such model predicts that the lineshape will be indepen- 
dent of vibrational transition, and that the R(J) and 
P(J — 1) transitions should be mirror images, which is 
not what is generally observed. Once vibrationally av- 
eraged potentials for different vibrational states become 
available, it will be possible to more critically test the 
lineshapes predicted by this model. At present, it must 
be admitted we do not have a quantitatively successful 
theory. 



C. Line shapes of rotationally resolved lines 

Direct evidence of inhomogeneous broadening is only 
available for a few molecules; nevertheless, there are rea- 
sons to believe that it is common to all molecules. One 
reason is that the ro-vibrational lines rarely fit to the 
Lorentzian shapes predicted by simple lifetime broaden- 
ing models. This is particularly evident for HON, for 
which the R(0) ( J = — > 1) transition has been studied 
in several vibrational states, as well as in the deuterated 
molecule (see Figure ||). For OCS, the Hnes in the P and 
R branch have tails that point towards the band origin 
and widths that grow with J. Possible explanations have 
been discussed by Grebenev et al. Available evidence 
suggests that one can expect inhomogeneous broaden- 
ing to contribute a few hundred MHz or more to the 
linewidth, equivalent to a lifetime of a few ns. Often, 
rotational and vibrational relaxation is slower than that, 
yet (with the notable exception of vibrationally excited 
HF [Q) still much faster than the 0.1-1 ms time scale 
of the experiment. 

This slow relaxation derives from the large mismatch 
between the frequencies of the modes of the molecule and 
those of the helium bath (the helium cannot respond on 
the time scale of the vibrational motion of the molecule). 
The maximum energy of an elementary excitation in bulk 
''He is about 20 K As a result, V-T (vibration to 

translation) and even most V-V (vibration to vibration) 
relaxations require simultaneous creation of many (10~ 
100) phonons/rotons in the helium. One expects the rate 
of such processes to be slow, decreasing exponentially 
with the number of phonons involved. 



1. Homogeneously broadened profiles 

Exceptions to the above argument are expected in 
cases of very low frequency, large amplitude modes, and 



in cases where there is a resonance, or near-resonance, 
of vibrational levels in a molecule and the optical ex- 
citation accesses an upper component of the resonance 
polyad. This can be compared to the behavior in the 
gas phase, where it is well established that anharmonic 
resonances allow for rapid (i.e, approaching the speed 
expected for rotational relaxation) V-V coUisional relax- 
ation inside the set of coupled levels. 

A related phenomenon is the case of molecules which 
in the gas phase have a density of states slightly below 
that required for IVR via high-order relaxation, such as 
propyne in the 2vi mode. In helium, the effective density 
of states is raised, as the helium can be excited with one 
or more low energy quanta, thus accommodating small 
energy mismatches. Molecules with a sufficiently high 
density of states to undergo statistical IVR in the gas 
phase, such as (CH3)3SiCCH are at most weakly 
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FIG. 5: R(0) line for DCN {vi = 0) and HCN (v^ = 0, 1,2), 
evidencing the importance of inhomogeneous broadening. 
These data show that vibrational excitation is not very im- 
portant in determining the broadening pattern, whereas the 
mass of the dopant is, suggesting a nontrivial dynamics inside 
the droplet. The data are adapted from Refs. |7^, |6^, 
respectively. The zero of the frequency was selected to best 
overlap the transitions. 
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affected by the helium. 

An interesting case occurs when, in the gas phase, ex- 
citation would cause the molecule to dissociate. While 
the helium forms a "cage" which may prevents dissocia- 
tion, it would still accommodate some of the excitation 
energy. As in the case above, one might still be left with 
an excited species, but as for determining the lifetime, 
the initially excited state would appear to have relaxed. 
So far this scenario has only been observed for van der 
Waals complexes, and the lifetimes are only moderately 
decreased relative to their gas phase value. All the values 
measured to date are reported in Table |l|. 

All the experimental evidence currently available indi- 
cates that the line profiles are determined by rotational 
relaxation only when the energy of one rotational quan- 
tum is sufficient to excite bulk modes of the droplet [|150[ . 
As the energy must then be near or above the roton min- 
imum of 8.58 K Q, only small molecules, with large ro- 
tational constants, display this kind of behavior. The ob- 
served lines are Lorentzian, with widths around 1 cm~^, 
reflecting relaxation times in the 1-10 ps range (see Ta- 
ble 



TABLE II: Vibrational lifetimes calculated from apparent ho- 
mogeneous linewidths in He nanodroplets (except Ar-HF, 
measured directly) . The entry under the header "Note" refers 
to the cause of vibrational relaxation; "tentative" indicates 
that vibrational relaxation is assumed as the dominant cause 
of broadening. 



Molecule Mode Lifetime/ns 



Note 



Ref. 



HCOOH 



di-HCOOH 

H"C"CH 
(HCN)2 
(HCN)2 
(HF)2 
(HF)2 
(HF)2 
Ar-HF 
SFe 
CH3CCH 



ui (a) 
'^1 (b) 
ui (c) 
1^1 
vs 

V3 
Vl 
V2 
V2 

1^2 + 1^5 

Vl 
V 



0.12 
0.09 
0.05 
0.15 
0.12 
« 1 

> 0.64 
0.15 

> 0.53 
7 X 10"^ 
3 X 10"^ 
1.5 X 10^ 

0.56 
0.16-1.5 



Fermi resonance 



Fermi resonance 
Fermi resonance 
Fermi resonance 

dissociation 
dissociation 
tentative 
tentative 

tentative 
tentative 



[147 

65 

65 
65 
65 

H 

28 
69 



TABLE HI: Rotational lifetimes measured in He nan- 
odroplets. The entry in the column "Note" refers to the 
droplet modes responsible for rotational relaxation. 



Molecule 



Levels Lifetime 
or band (ns) 



Note 



Ref. 



HCN/DCN J = ~ 10 surface? 
HCCCN J = 3, 4 2-20 

J = 2, 4 10-100 
HF i/i 12 X 10~^ bulk 

(HF)2 VI 2.7 X 10"^ " 



^1 

68] 

ii 

65] 



2. Indirect measurements: Saturation, double resonance 

For large molecules, the dissipation channel into bulk 
helium modes is closed for thermally populated rota- 
tional states, and relaxation must occur either by R-T 
(rotation-translation) transfer (for which Franck-Condon 
factors are unfavorable |^) or by coupling to the softer 
surface modes of the droplet. Because long-range forces 
confine the molecule to spend most of its time away from 
the surface, the latter coupling is weak, and relaxation 
times are corres pon dingly extended into the nanosecond 
range (see Table III). In fact, the coupling to the surface 
modes has been estimated to be so weak that it is 
an outstanding question how the molecules can come to 
equilibrium with such a grainy heat bath. Lines are then 
inhomogeneously broadened, but homogeneous widths 
can still be estimated from saturation and double res- 
onance measurements. 

The first such measurement, the MW-MW experiment 
on HCCCN [Q, faced a set of apparently contradictory 
findings: the single resonance saturation data demon- 
strated that the strength of rotational lines was propor- 
tional to the square root of the applied power over a wide 
range, the classic signature of saturation of a line dom- 
inated by inhomogeneous broadening As the only 
previously considered source of inhomogeneous broaden- 
ing was the spread of droplet sizes, a standard analysis 
based on the assumption of "static" broadening predicted 
that double resonance measurements would directly re- 
veal the homogeneous component of the line via hole- 
burning with the the pump field, i.e., a Bennet hole and 
hill in the absorption of the probe [7^, with a width of 
no more than 150 MHz. It was instead found that the 
linewidths of the transitions affected by the pump beam 
were comparable to those measured in single resonance 
(1 GHz). The contradiction was solved by ascribing in- 
homogeneous broadening to the spectral diffusion of the 
dopant over a set of quantum states (most likely different 
translational states within the droplet; see Section III E| ) . 
Initially it was assumed that such diffusion occurred on 
a time scale faster than rotational population relaxation. 
Upon further analysis, it was realized that this would 
have corresponded to a much higher saturation power 
than was actually observed, and it was proposed that 
the two phenomena are occurring on the same time scale 
(tens of ns) |6^. Despite the use of a simplified model, 
and the tentative value of the parameters obtained from 
it, the conclusion that spectral diffusion and rotational 
relaxation both occur because of the motion of the dopant 
inside the droplet appears to be reasonable; indeed the 
time scale for these is in the ns range. 

MW-IR experiments performed on HCCCN Q and 
OCS 1 80 confirmed the above picture: double resonance 
lines span the whole inhomogeneous width. A repro- 
ducible substructure of 5 peaks, about 50 MHz wide is 
observed for the J — 2 i rotational transition of OCS, 
suggestive of a 2 J -I- 1 splitting, but no pattern is re- 
solved for other lines, thus making a definite assignment 
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not possible. The main finding, common to both experi- 
ments, and therefore probably of general nature, is that 
the populations of all the rotational levels are altered by 
the MW "pump" (Figure |). 

In gas phase DR experiments, such as the classic 
IR-MW studies by Oka one can often observe 

collisionally-induced four level DR signals, but they are 
typically considerably weaker than the coUisionless fea- 
tures. In helium, the relaxation induced signals have 
been found to be even stronger than the "direct" sig- 
nals in some cases. There also was a strong dependence 
of the DR signals on the source conditions, which sug- 
gests that the relaxation dynamics are strongly depen- 
dent upon droplet size, which affects the density of states 
of the thermal bath that the molecule is coupled to. For 
HCCCN, the DR spectra dramatically simplify for the 
largest droplets, which has suggested rotational relax- 
ation is in the "strong collision" limit for these droplets, 
i.e., that the relative populations produced by relaxation 
are independent of the initially populated state. 

A curious feature observed for both HCCCN and OCS 
is that there appears to be a lack of conservation of pop- 
ulation in the DR experiment. For DR signals that arise 
from population transfer, the integral of the DR signal 
over the rotation structure of the IR band should be zero 




2061.2 2061.4 2061.6 2061.8 2062.0 2062.2 



Wavenumber, [cm"^] 

FIG. 6: Double resonance spectra of OCS in He droplets. 
The levels coupled by the MW "pump" and the frequency of 
the transition are indicated in the upper left corner of each 
panel. In contrast to what is commonly observed in gas phase 
experiments, note how all the levels are affected by the MW 
"pump" , and how the total signal does not sum to zero. From 
Ref. g. 



since the integral cross section for IR absorption is inde- 
pendent of the lower state. This apparent lack of pop- 
ulation conservation must reflect some "artifact" in the 
experiments, but the ones considered do not appear con- 
sistent with the observations; in particular, the sign of 
the total signal is opposite for the two molecules. This 
effect is currently unexplained and may ultimately reveal 
some interesting novel dynamics for these systems. 



D. Molecular rotations 

1. Changes in rotational constants: How is the helium 
moving? 

An important aspect of the rotational structure of sol- 
vated molecules is that in most cases the effective rota- 
tional constants are substantially reduced from their gas 
phase values. For molecules with rotational constants 
much larger than 1 cm~^, the change on going into he- 
lium is small, while for heavier molecules, the rotational 
constant is typically reduced by a factor of 3-4 upon he- 
lium solvation. Clearly, when the molecules rotate in he- 
lium, at least some of the helium must also move, creating 
kinetic energy and contributing to the effective moment 
of inertia of the molecular rotation. Several models have 
been put forward to describe the nature of the helium 
motion. 

The first model that was proposed Q assumed 
that a specific number of helium atoms were rotating 
rigidly with the molecule, directly contributing their mo- 
ment of inertia to that of the molecule. This model ap- 
peared to work for both SFg and OCS, in that chemically 
reasonable positions and numbers of rigidly attached he- 
lium atoms were consistent with the experimentally ob- 
served changes in moment of inertia (A/). For OCS, the 
authors called this the "donut model" since it assumes 
that a ring of 6 He atoms is localized in the potential 
minimum locus around the OCS molecule. This very 
simple model has two problems. First, when one tries 
to extend it to lighter molecules such as NH3 or HCN, 
even a single rigidly attached He atom would increase the 
moment of inertia by nearly an order of magnitude more 
than what is observed experimentally. Second, more re- 
cent results for the OCS-H2 complex in helium are 
clearly inconsistent with the "donut model" for OCS. In 
this work, it was found by isotopic substitution that a 
single H2, HD, or D2 molecule will replace one of the he- 
lium atoms in the "donut" . Yet, extrapolation to zero 
mass of the isotopomer gives a moment of inertia of the 
"donut with a bit taken out" considerably larger than for 
the OCS in a pure Helium droplet, i.e. with a complete 
"donut" . 

A more sophisticated "two fluid" model for the motion 
of the helium was introduced by Grebenev et al. 
They proposed that the density of helium around the 
molecule be partitioned into spatially dependent normal 
and superfluid fractions, with the normal fraction having 
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large values only in the first solvation layer. In analogy 
with the classical Andronikashvili experiment (where the 
normal fluid fraction of bulk helium follows the rotation 
of a stack of pla tes, contributing to their effective mo- 
ment of inertia |100| ), these authors propose that the 
normal fluid fraction rotates rigidly with the molecule, 
contributing the classical moment of inertia of the normal 
fluid mass density. This model was incomplete, however, 
in that the authors provided no definition of the spatially 
dependent normal fluid density that appears in the ex- 
pression for the increase in moment of inertia. 

Kwon and Whaley proposed an ansatz for a local su- 
perfluid density |101|, based upon an extension of Ceper- 
ley's global estimator for the bulk superfluid density 
from path integral Monte Carlo (PIMC) calculations 
They took the difference in total density and superfli 
density to be a nonsuperfluid density that rotates rigidly 
with the molecule, assuming the molecule-Hc interaction 
is sufficiently anisotropic. The term "nonsuperfluid den- 
sity" is used instead of "normal fluid density" because 
this density is not related to the normal fluid density 
produced by a "gas" of elementary excitations, as in bulk 
helium experiments such as the Andronikashvili experi- 
ment. The inertial response of the helium to rotation is 
a second rank tensor, involving both the axis of rotation 
and the direction of the response, which is measured in 
PIMC by a projected area of closed Feynman paths [ p9| . 
The Kwon- Whaley estimator produces a scalar superfluid 
density, and thus has the wrong symmetry to describe 
the helium inertial response [102, 103|. This two fluid 



uid 



model has, however, made predictions in excellent agree- 
ment with experiment for the cases of SFg and OCS in 
helium. Draeger and Ceperley have proposed a new local 
superfluid density estimator that has the correct tensor 
properties and gives, when spatially ave raged , the correct 



value for the global superfluid fraction 103 



An alternative, quantum hydrodynamic model has also 
been developed for the helium contribution to the mo- 
ment of inertia of molecules in helium. It is well known 
that linear motion of a rigid body through an ideal (avis- 
cous and irrotational flow) fluid generates motion in the 
fluid that contributes to the effective mass for transla- 
tion. For example, for a sphere, the effective mass in- 
crease is exactly one half of the mass of the displaced 
liquid |104|. Similarly, rotation of an ellipsoid generates 
helium kinetic energy proportional to the square of the 
angular velocity, thus contributing to the effective mo- 
ment of inertia. The possible inclusion of such a "su- 
perfluid hydrodynamic" term was already recognized by 
Grebenev et al. [^. However, attempts to estimate this 
term, using classical expressions for the moment of in- 
ertia of an ellipsoid in a uniform fluid of the density of 
bulk helium, gave results which were only a small frac- 
tion of the experimental values 



91 . On the other 



hand, as the density around a molecule is highly nonuni- 
form, Callegari et al. |105| have developed a quantum hy- 



inertia in good agreement with experiments for a number 
of molecules. 

If one adiabatically separates the helium motion from 
the rotation of the molecule, then one expects the helium 
density in the rotating frame of the molecule to be con- 
stant and equal to that around a static molecule |10E, 
106|. This is called the adiabatic following approxima- 



tion and, on physical grounds, should be accurate when 
the minimum rotational spacing, 25, is well below the 
energy of the roton, which characterizes the helium den- 
sity excitation energy on the atomic level (in the two fluid 
model, Kwon et al. ha ve introd uced a distinct definition 
of adiabatic following |101, 107|). If one further assumes 
that the fluid velocity field generated by this time de- 
pendent helium density in the laboratory frame is irrota- 
tional (which Kelv in proved to be the lowest kinetic en- 
ergy solution [108|), then one obtains from the equation 
of continuity a second order differential equation for the 
velocity potential (the instantaneous laboratory frame 
velocity is the negative gradient of this velocity poten- 
tial). The helium contribution to the moment of inertia 
is then calculated directly as an integral over the velocity 
potential. In applying this method to molecules of cylin- 
drical symmetry, the helium density was calculated us- 
ing a density functional (DF) approach |109, IIC] which, 
instead of searching for the ground state of the helium 
many-body problem, finds the solution of a one-particle 
problem with a nonlocal self interaction, whose density 
distribution, if the correct functional was used, would ex- 
actly match the density of the many-body wavefunction. 




0.01 



0.1 



10 



drodynamic model that properly takes into account the 
anisotropy of the density, and have obtained moments of 



FIG. 7; Plot of the rotational constants (gas phase versus 
helium droplet) reported in Table ^ The two lines are the 
functions y = x and y — x/3 respectively. Note the differ- 
ent trend above and below Bgas ~ 1 cm~^, indicative of a 
threshold for adiabatic following. 
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It takes orders of magnitude less computational effort 
to solve for the helium density by DF methods than by 
Monte Carlo methods, and the results have been found 
to be in good agreement for cases where the molecule- 
helium interaction potentials are not too deep [111|. 



An alternative derivation [112| of the quantum hy- 



drodynamic approach is to write the many-body helium 
wavefunction in a rotating "external" molecular potential 
as that of the instantaneous static ground state wave- 
function times a single particle phase function. Varia- 
tional optimization of the energy with respect to this 
phase function gives the same hydrodynamic equation as 
for an ideal classical fluid with the adiabatic following ap- 
proximation. Within a Hartree approximation for the he- 
lium wavefunction (which DF theory implicitly invokes), 
the hydrodynamic treatment gives the exact wavefunc- 
tion correction to first order in angular velocity and thus 
the energy to second order. Time reversal symmetry im- 
plies that changes in the helium density can only depend 
on even powers of the angular velocity, and will therefore 
contribute only to the centrifugal distortion constant, not 
to the rotational constant. It has been known since the 



dawn of quantum mechanics |113, 114| that quantum 
mechanics can be formulated exactly as a hydrodynamic 
problem in configuration space. The one body approxi- 
mate hydrodynamic approach used for these calculations 
is equivalent to many calculations that have been done 
in bulk heli um, including the structure of the core of a 
vortex 1 115 1, which has a helium density structure on a 
length scale smaller than that induced by the solvation 
of most molecules. 

The quantum hydrodynamic model has been applied to 
a number molecules, including HCN, HCCH, HCCCII3, 
OCS, HCCCN, and HCN dimer IT05|. For the last four, 



for which the adiabatic following approximation is ex- 
pected to hold, the predicted increase in moment of in- 
ertia was between 103-130 % of the observed values, as 
good as could be expected given the uncertainties in the 
helium solvation densities. For HCN and HCCH, which 
have small moments inertia, the hydrodynamic calcu- 
lations significantly overestimate the observed increase 
which is taken as evidence of the breakdown of the adi- 



abatic following (see below). Kwon et al. [107| have re- 
ported that the hydrodynamic model predicts for SFg an 
increase in the moment of inertia of only 6(9)% of the ex- 
perimentally observed value when their calculated super- 
fluid (total) helium solvation density was used. Similar 
hydrodynamic calc ulations have been done by Lehmann 
and Callegari |116] using the lowest four tensor compo- 
nents in the expansion of the helium density in octahedral 
harmonics. These authors find that when the total he- 
lium density is used, the hydrodynamic contribution to 
the moment of inertia was 55% of the observed value, 
and this value is likely to increase if the full anisotropy 
of the helium density is included. The reason for the 
discrepancy in the two reported values |107, |116[ for the 



Despite their apparent differences, the "two fiuid" and 
the hydrodynamic approaches to the superfluid helium 
rotation problem are related. They are both based on 
a linear response of the helium, and thus not appli- 
cable when the adiabatic separation of molecular rota- 
tion and helium motion breaks down. Furthermore, the 
"nonsuperfiuid" helium accounts for much of the angu- 
lar anisotropy of the helium density, which is the source 
term in the quantum hydrodynamic treatment. 



hydrodynamic prediction for SFg has not yet been re- 
solved. 



Clearly, it would be desirable to develop a theoreti- 
cal method to directly calculate the rotational constant 
of molecules in superfluid helium without making the dy- 
namical approximations of either the "two fluid" or quan- 
tum hydrodynamic approaches. It would also be desir- 
able to directly calculate relaxation rates for excited rota- 
tional states as well as their energies. Unfortunately, this 
is still largely an unsolved problem. Blume et al. |117l 
have used an excitation energy estimator based upon 
the rate of decay of an approximate excited state wave- 
function under imaginary time propagation (POITSE 
method). Unfortunately, only studies of very small clus- 
ters have been published to date with this method, and 
thus its relevance to rotational dynamics in helium nan- 
odroplets is not established. Further, the initial state had 
the same rotor based nodal surfaces and thus the accu- 
racy of this method may also depend upon an approxi- 
mate adiabatic separation of the rotational motion from 
that of the helium. If this approximation breaks down, 
the initial state will contain sizable contributions from 
several excited states and the imaginary time propaga- 
tion will contain sizable c ontri butions from many expo- 
nential decays. Lee et al. [ [106| have calculated the rota- 
tional constant of SFg-Heyv {N — 1-20) by a fixed-node 
diffusion Monte Carlo method. In these calculations, the 
helium and the molecule moved in the same coordinate 
system, a condition that is particularly easy for a spher- 
ical top molecule, but harder in general. Calculations 
for the rotational spacings of SFg were in quantitative 
agreement with experiment for A'^ > 8, as in the simple 
model proposed earlier by Hartmann et al. p3| The fixed 
node approximation used, based upon the nodes of the 
rotational wavefunction of the isolated molecule, was jus- 
tified by precisely the same type of adiabatic separation 
of time scales that is used for adiabatic following in the 
quantum hydrod ynam ic approach [106|. Recent work on 
a model system | 118 | has found that the hydrodynamic 
and this molecular rotor based fixed node approximation 
arc closely related, and break down at about the same 
point when the rotational constant of the molecule is in- 
creased. To date, the diffusion Monte Carlo calculations 
have been published only for SFg, so the accuracy of the 
method in general is hard to judge, though it is likely the 
most accurate of the presently available methods. 
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2. Breakdown of helium adiabatic following 



As it was mentioned in the beginning of Section III D 1 



"light" molecules with large gas phase B values appear to 
have little helium contributions to their moments of iner- 
tia. It is certainly the case that these molecules, by their 
nature, have relatively weakly anisotropic interaction po- 
tentials with the helium, and thus a small increased mo- 
ment of inertia is expected. However, because they also 
have very small molecular moments of inertia, even a 
small helium contribution would be expected to lead to 
a substantial fractional change in the value of B. This 
suggests that something else is going on. The models dis- 
cussed above invoke "adiabatic following" of the helium 
with the molecular rotation, though the precise meaning 
of that phrase is different for the quantum hydrodynamic 
and the "two fluid" models. It is expected that this adi- 
abatic separation of the time scales of molecular rotation 
and helium motion will break down when the rotational 
spacing becomes comparable to the energy required to 
create a helium excitation with a wavelength on the order 
of the solvation structure around the molecule, and this 
will be similar to the energy of the helium roton, 8.58 K 
|p8| . This suggests a rotational constant of w 4.3 K as the 
point where substantial breakdown of the adiabatic fol- 
lowing would be expected (see Figure |7|). Lee et al. 106 



studied adiabatic following breakdown theoretically by 
changing the rotational constant of SFg in their fixed 
node DMC calculations of small SFg-HeN clusters. They 
found that the angular anisotropy of the helium in the 
molecular frame and the helium contribution to the ef- 
fective moment of inertia both decreased substantially 
when the SFg rotational constant used in the calcula- 
tion was increased by a fa ctor of ten, from 0.1 Icm"^. 
HCN {B = 2.13 K ||ll|l) is one of the few molecules 
that have been studied which is near this threshold. The 
breakdown of adiabatic following for HCN has been es- 
tablished experimentally, by observation of the change in 
the helium contribution to the effective moment of iner- 
tia by isotopic substitution of HCN to DCN. Compared 
to DCN, it was found that the somewhat lighter and thus 
23% faster HCN produced a 9% decrease in the helium 
contribution to the moment of inertia |Q. It would be 
highly desirable to study other molecular systems in this 
important intermediate region. 

Density-functional computations can reproduce the 
breakdown of adiabatic following to some degree by ex- 
plicitly moving the density computation from the inertial 
lab frame to the frame rotating with the dopant molecule. 
In terms of the helium and total angular momentum op- 
erators L and J, respectively, the rotating-frame rota- 
tional energy of the dopant molecule is B{J — L)^ = 



BJ^+BL'^-2BJ-L |120| . In the rotational ground state, 
(J^) = 0, and (L) = because of time reversal symme- 
try. Therefore, the only change in going to the rotating 
frame is the appearance of a term BL^ in the helium den- 
sity functional. For large values of B, this term increases 



the angular smoothness of the ground state helium den- 
sity, which decreases the helium moment of inertia in the 
hydrodynamic theory; this provides a scale to gauge the 
breakdown of adiabatic following. Work is in progress in 
our group to quantitatively assess this approach. 



3. Giant centrifugal distortion constants 

Another interesting generic feature that has been ob- 
served in rotational spectra in helium droplets is a dra- 
matic increase in the centrifugal distortion constants re- 
quired to fit the spectra. Values of D/B ^ 10"'* have 
been observed for heavy molecules, three to four orders 
of magnitude larger than for the same molecules in the 
gas phase. It has been suggested that these large distor- 
tion constants reflect weakly bound helium atoms that 
are easily displaced by the centrifugal potential, but a 
careful attempt to estimate this effect, for OCS, gave 
a predicted D value 50 times smaller than is observed 
experimentally js^. If one uses the standard scaling 
D ~ AB^/u!^, one must invoke an effective vibrational 
wavenumber on the order of 1 cm^^ in order to repro- 
duce the observed D values, much smaller than what is 
physically reasonable considering the force constants of 
the helium-molecule potentials. Given the above discus- 
sion of the breakdown of adiabatic following with increas- 
ing rotation rate, one can argue that the effective D value 
should be negative, the opposite sign from what has been 
observed! This is because we would expect the effective 
B value to increase, going towards the gas phase value, 
as the rotational velocity increases, making it harder for 
the helium to follow. 

At present, there is only one calculation that has re- 
produced t he c orrect size of the distortion constant: a 
toy model |llS| of a rigid, planar ring of N He atoms in- 
teracting with a rotating molecule. What this calculation 
clearly demonstrated is that opposite behaviors are found 
depending on how the rotational velocity is increased: by 
decreasing the moment of inertia of the molecular rotor 
or by spinning up a fixed rotor to higher J values. In 
the former case, the helium follows more poorly, consis- 
tent with the above discussion. In the latter case, the 
He anisotropy and following increase until a resonance 
between the molecular and helium ring rotations occurs. 
These results are consistent with general conclusions pre- 
viously drawn by Leggett in his analysis of the "rotating 
bucket" experiment for bulk liquid helium [121|. 



4- Vibrational dependence of rotational constants 

As in the gas phase, the effective rotational constants 
change upon vibrational excitation. The size of the ob- 
served vibration-rotation coupling is much larger in he- 
lium than that expected based upon the gas phase change 
in B. For example, for HCCCN the total increase in mo- 
ment of inertia upon vibrational excitation changes from 
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0.35 uA^ in the gas phase to 9.2 uA^ in hehuni, and for 
CH3CCH it changes from 0.32 uA^ to 8.5 uA^ Q. This 
suggests that most of the effect is due to a change in the 
hehum contribution to the effective moment of inertia. 
As expected, the excited vibrational states have increased 
hehum contributions, which suggests increased density 
and/or anisotropy around the vibrationaUy excited state. 
As pointed out above, the absence of phonon wings to the 
spectra indicate that these changes must be smah. There 
has yet to be pubhshed a quantitative attempt to cal- 
culate this change in solvation structure, which requires 
knowledge of the vibrational dependence of the molecule- 
He interaction potential, which are not yet generally 
available. Recently, there has been important progress 
made in the calculation of the intermolecular-mode de- 
pendence of the interaction potentials |122{ , which would 
allow such an estimate to be made. 



E. Shift in band origins 

An extremely attractive feature of ro- vibrational spec- 
troscopy in helium is that solvent shifts of the vibrational 
origins are generally quite small, typically less than 0.1%. 
This can be contrasted with spectroscopy in the tradi- 
tional matrix environments of Ne or Ar, for which the 
shifts are usually much larger. As an example, for HF 
the m atrix shifts are 0.067%, 0.24% and 1.1%, respec- 
The small shifts in helium are due in large 



tively [123 



part to its low density and polarizability (even when com- 
pared to other rare gases). However, existing evidence 
suggests that there is often a fortuitous cancellation of 
effects leading to smaller shifts than one could otherwise 
expect. 

In considering the shift in band origin, one must dis- 
tinguish between the vertical and adiabatic excitation en- 
ergies. The former refers to a transition within a frozen 
helium solvation density, while the latter considers the 
optimized helium solvation density in each vibrational 
state. As discussed above, the observed vibrational de- 
pendence of the effective rotational constants is largely 
a result of such a change in helium solvation. The peak 
of the zero phonon line should provide an experimen- 
tal estimate for the adiabatic excitation energy while the 
vertical excitation will be at the 'center of gravity' of the 
entire transition, including the phonon wing. The differ- 
ence between vertical and adiabatic excitation energies 
reflects the helium reorganizational energy during vibra- 
tional excitation. Given the lack of observed "phonon 
side bands" in vibrational spectra, one can anticipate 
that this reorganizational energy is small, but not nec- 
essarily small compared to the typically observed helium 
induced vibrational frequency shift. If a ~ 10 cm^^ wide 
phonon wing contains 50% of the integrated intensity of 
a transition, it would cause a ~ 5cm~^ shift between 
the adiabatic and vertical excitation wavenumber, but 
this phonon wing would have a peak absorption strength 
less than 1% of that of the zero phonon line (assuming 



a linewidth of 1-3 GHz), and thus could be easily lost in 
any low frequency noise in the spectrum. The vertical 
excitation energy, in contrast, is much easier to estimate 
theoretically, since it does not require computation of the 
change in helium-helium interactions. As an example, 
Blume et al. [|5^ have calculated the vertical excitation 
wavenumber of HF in helium by averaging the potential 
difference for He interacting with v=0 and 1 states of HF 
over the ground state wavefunction of the system. 

The change in the helium- molecule interaction poten- 
tial upon vibrational excitation reflects several differ- 
ent physical components. Based upon the natural sep- 
aration in time scales between intramolecular and he- 
lium solvation motions, an adiabatic separation is ap- 
propriate and one can consider a vibrationaUy aver- 
aged molecule- heliu m interaction potential for each vi- 
brational state |l22| . If we consider excitation of a nor- 
mal mode fc in a molecule with dimensionless normal 
coordinates q^, then in the harmonic approximation the 
He-molecule potential in the excited state of the normal 
vibration can be written as 



dV 



dql 



where V is the full interaction potential, including in- 
tramolecular coordinates, and the derivatives are evalu- 
ated at the equilibrium structure of the isolated molecule. 
The sum is restricted to totally symmetric vibrational 
modes, and (qs) is proportional to the cubic anhar- 
monic coupling (jfskk between the excited mode k and 
the symmetric mode s. Given the small changes in vi- 
brational wavenumber upon solvation by helium, the har- 
monic approximation is expected to be adequate for rigid 
molecules. 

The long range interaction between He and the 
molecular dopant is dominated by the isotropic and 
P2{cos6)R~^ terms. One term contributing to the vi- 
brational dependence of this term is the polarization sta- 
bilization of the dipole induced by vib ratio nal motion, 
even in nonpolar molecules such as SFg |124|. This term 
leads to a shift proportional to the IR strength of the 
vibrational mode. Integrating this term over an estimate 
for the helium solvation density gave a shift for SFg 
of about one half the observed value of —1.6 cm~^ p4| , 
(such an estimate is for the vertical, not adiabatic shift). 
A sizable fraction of the shift comes from the flrst sol- 
vation shell, for which the use of the long range 
terms is not expected to be complete, though there may 
be a signiflcant cancellation of errors between the ne- 
glect o f hig her order terms and damping of the dispersion 
series [ 126 |. For molecules with a permanent dipole mo- 
ment, one also needs to consider the changes in the vibra- 
tionaUy averaged dipole moment. While these changes 
are usually smaller than the transition moments, at least 
for reasonably strong IR transitions, one gets an "inter- 
ference" with the permanent dipole moment (since the in- 
duction energy depends on (/i^)), which can magnify the 
effects of a small absolute change in the dipole moment 
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upon vibrational excitation. For excitation of the over- 
tone of the C-H stretching mode in HCN |^ , the tran- 
sition dipole and change in vibrationally averaged dipole 
moment lead to estimated shifts of 0.14 and 1.67 cm~^, 
respectively. The long range interaction between He and 
molecular dopants is typically dominated by dispersion 
interactions, not the induction terms. One can estimate 
the change in the induction terms from changes in the po- 
larizability of the molecule upon vibrational excitation. 
For HCN, this gives an estimate of 4.98 cm~^ for the 
shift, substantially larger than the observed value. 

All of the terms discussed above are long range and 
lead to red shifts. They certainly are expected to dom- 
inate the droplet size dependence of the shift. Putting 
the solute in the center of the droplet and integrating 
the missing (compared to the bulk) interaction, one finds 
that the shift from bulk or infinite droplet size limit 
will be dominated by the terms and will scale as 
N~^, where N is the number of helium atoms in the 
droplet |3l|. Thus, a red shift towards the limiting 
bulk value is expected for large droplets, independent 
of whether the overall shift is red or blue. If one as- 
sumes that the droplet size distribution follows the log- 
normal distribution, then the inhomogeneous distribu- 
tion of droplet sizes will also contribute a log-normal dis- 
tribution of shifts away from the bulk value. This al- 
lows the N dependence of the width and shift for smaller 
droplets (where the inhomogeneous cluster size distribu- 
tion contributes most to the spectral shape) to be used to 
estimate the inhomogeneous size dispersion contribution 
to the lineshape of larger droplets. 

Counteracting the above discussed increases in long 
range He attraction upon vibrational excitation (and 
higher order multipole contributions that have not yet 
been estimated) is an expected blue shift from the 
"crowding" interaction. Such an effect is due to the fact 
that the potential for a molecular oscillator will be in part 
stiffened due to the interaction with the first solvent shell, 
where repulsive forces are important. In a highly simpli- 
fied harmonic picture of a linear X-H-He unit, the X-H 
stretch normal mode force constant will be increased by 
the H-He force constant. This leads to a blue shift of 
the vibrational frequencies that should be most impor- 
tant for vibrational motion dominated by exposed atoms 
and/or involving low frequency, large amplitude motions. 
Exa mple s of the later include the torsional mode of gly- 
oxal 1 127] and the butterfly mode of pentacene j7|], both 
in the first excited electronic states. At fixed helium sol- 
vation structure, one would expect that the blue shift of 
a chemically unique vibrational motion, such as the C-H 
stretching mode of terminal acetylene compounds, would 
be approximately conserved. However, helium density 
functional theory calculations made by our group have 
shown that helium solvation density around this terminal 
C-H bond is a strong function of the rest of the molecule, 
which thus far has prevented the development of a quan- 
titatively predictive theory of the expected blue shift. 

The spectrum of HCN dimer provides some evidence 



for the importance of the blue shift induced by the local 
structure. For t he hydrogen bound C-H stretch, Nauta 
and Miller |12S| have observed a larger gas-to-helium- 
droplet redshift compared to the same quantity for the 
HCN monomer. They interpreted this as arising from an 
increasingly linear and thus strengthened hydrogen bond 
when HCN dimer is solvated in helium. However, this ef- 
fect can be also interpreted as arising from a substantial 
reduction of the helium induced blue shift upon inter- 
molecular complex formation since the hydrogen bonded 
H atom no longer pushes directly on the helium solvent. 

At present, we lack an effective theory to use the ob- 
served solvent induced shifts to learn about the solvation 
structure around the chromophore molecule or to reliably 
predict the observed shifts, even their sign, relative to the 
gas phase. However, since the observed shifts appear to 
be small, we can safely neglect this complication in most 
cases, in particular when using the predictions of quan- 
tum chemistry to help assign novel spectra observed in 
liquid helium. 



IV. SPECIAL TOPICS 

A. Synthesis of nonequilibrium structures 

A very promising application of HENDI spectroscopy 
which has been recently and dramatically demonstrated 
by Nauta and Miller ^ [gJ] relates to the possibil- 
ity of using the very cold droplet environment to syn- 
thesize highly non-equilibrium species that could not be 
prepared in any other known way. In analogy to the 
preparation of van der Waals clusters of spin polarized 
alkali atoms demonstrated in Princeton in 199 6 (se e the 
following review by Stienkemeier and Vilesov |12£| and 



Refs. 130, 131, |132D , the essential ingredient in these 
syntheses is, in addition to the low temperature of the 
droplets, the rapid cooling action that the superfluid he- 
lium exerts on the molecular partners at the time of 
cluster formation. Consider the example of the prepa- 
ration of linear chains of HCN |^: such chains with 
more than three monomer units have not been observed 
even in the cold environment of a supersonic jet expan- 
sion, where cyclic or other nonpolar structures are in- 
stead prepared. When an HCN molecule collides with a 
cluster that already hosts another molecule in its interior, 
the long range dipole-dipole forces align the two dipoles 
at distances where no other force is acting between the 
two molecules. The rapid rotational relaxation necessary 
for such an alignment is probably provided during mul- 
tiple collisions with the droplet surface, where coupling 
to ripplons is enhanced. After the two dipoles have been 
aligned, the molecules are drawn together by their mu- 
tual attraction; in contrast to a gas phase coagulation, 
however, their potential energy is not completely trans- 
formed into kinetic energy because of the steady cooling 
action of the "solvent" . Consider that a potential energy 
of 200 K (i.e., 2.8 x 10~^^ J) would accelerate a molecule 
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of mass 27 u to a velocity of 350 m/s , which is weh above 
the critical velocity of 60 m/s |133] above which a for- 
eign body will rapidly dissipate energy in superfluid he- 
lium. This dissipative cooling during complcxation en- 
sures that, at the contrary of what happens in the gas 
phase, the freshly formed complex is born cold, without 
the energy needed to surmount any barrier to isomeriza- 
tion toward the global energy minimum. It is remarkable 
that this process is repeated until the linear chain spans 
the whole cluster diameter, as the data of Nauta and 
Miller show. Observed spectra for CH3OH and CH3CN 
clusters [p9|, for H2O clusters [p4| and for Ar„+HF clus- 



ters (described in their article in this issue 1 134 ) have 
been interpreted as arising from the isomers which are 
formed with the smallest amount of nuclear rearrange- 
ment. 

So far this self-assembling technique has been applied 
using stable molecules. Efforts are in progress in our 
laboratory to use CN radicals for preparing isocyanogen- 
like chains or, if the presence of a recombination barrier 
is detected, at least van der Waals clusters of CN radi- 
cals. The potential applications of this technique for the 
synthesis of exotic or high energy density materials are 
quite evident. 

A different, but equally useful, non-equilibrium syn- 
thesis has been demonstrated by Nauta and Miller in an 
experiment [|67| in which a cluster of reactive metal atoms 
(Mg) has been attached to an IR active molecule (HCN) , 
and the rotational spectrum of the resulting complex has 
been resolved (Figure ||). While clusters of metal atoms 
are bound by van der Waals or metallic forces accord- 
ing to their size [135, 136, 137, 138|, they are in both 



cases very difficult to study because their UV spectra are 
made of broad lines and their IR spectra, when present, 
appear at rather low frequencies, where lasers are not 
yet available. Probing the IR spectrum of a molecule at- 
tached to the atomic cluster yields valuable information 
on the latter via an analysis of the rotationally resolved 
spectrum or via the spectral shifts of the molecule's vi- 
brations. While the information contained in rotationally 
resolved spectra is at present clouded by our imperfect 
knowledge of the helium contributions (see Section [11 D| ), 
the type of spectrum (i.e., the symmetries of complexes) 
and shifts are already now sources of reliable information. 
We anticipate that this type of experiments will provide 
a wealth of information on metal-ligand interactions and 
on the nature of bonding in chemisorbed systems. 

Most large organic molecule, including those of bio- 
logical relevance, have a number of isomers populated 
at room temperature, particularly arising from rotation 
about single bonds (conformers) . When these molecules 
are cooled to low temperature, either in MI or super- 
sonic jet spectroscopies, isomers separated by high barri- 
ers do not interconvert, while those separated from lower- 
lying minima by low barriers will cool out. Multiple 
isomers have been observed by HENDI spectroscopy as 
well. As an example, see the study of glycine by Huisken 
et al. Byl, where the same three conformers were ob- 



served in He nanodroplets and in Ne and Ar matrices. 
In general, however, the different cooling rates should re- 
sult in different distributions of conformers in the differ- 
ent spectroscopies. Spectra of tryptophan and tyrosine 
in helium have been interpreted in terms of suppression 
of certain conformers in helium compared to those ob- 
served in seeded beam expansions [139 . The more rapid 



quenching, lower final temperatures, and shorter inter- 
action times of HENDI suggest that the conformer dis- 
tribution observed by this method should more closely 
resemble that of the initial ensemble. The much smaller 
spectral shifts in helium should make conformer assign- 
ments based upon quantum chemistry calculations sig- 
nificantly more reliable than for MI. The use of pendular 
spectroscopy in helium, which allows measurement of the 
permanent dipole moments and the angle between that 
and the vibrational transition moment, should further 
greatly enhance reliability of assignments. 



B. Large-amplitude molecular motions in liquid 
helium 



One of the great differences between electronic and 
ro-vibrational spectroscopy in He droplets is the ex- 
tent to which helium is displaced upon excitation of the 
molecule. Helium displacement arises from the repul- 
sive interaction (which can be viewed in terms of the 
Pauli exclusion principle) between the electron clouds of 
He atoms and that of the molecule. Electronic excita- 
tion results in general in an expansion of the electron 
cloud, with an energetic cost that is responsible for the 
experimentally observed blue shifts of electronic spectra 
of atomic and molecular species immersed in liquid He. 



Experiment 
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FIG. 8: Rotationally resolved spectrum of the van der Waals 
complex Mg3-HCN. The structural information extracted 
from this spectrum allows one to deduce a contraction of the 
HCN-metal-cluster bond, when going from Mg2 to Mga. This 
in turn indicates the presence of strong nonadditivity effects. 
From Ref. |^ 
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Conversely, the volume and shape of the electron cloud 
are not significantly changed upon ro-vibrational excita- 
tion. The observed band shifts are small, mostly towards 
the red, and are determined by a delicate balance be- 
tween the attractive and repulsive force s in the ground 
and excited state (see Section HIE ). As a first ap- 
proximation, one can say that molecular stretch vibra- 
tions, with few exceptions the only modes studied so far, 
are too stiff to be significantly affected by the presence 
of the helium. 

This statement does not hold true for large-amplitude 
modes which are, somewhat by definition, softer (the 
FWHM of the 1-D harmonic oscillator is proportional 
to (too;)^^/^ or {km)~^^^). They are particularly rele- 
vant not only because they displace more helium (sup- 
posedly coupling better to He bulk modes), but also be- 
cause they are usually bending modes. Tunneling barri- 
ers commonly occur for such modes, and are expected to 
be increased by the presence of the helium. For the above 
reasons, these mode are a useful probe of the molecule- 
helium interaction. 

The prototypical transition is the 1/2 umbrella mo- 
tion of NH3 [ |60| : the hindrance induced by the helium 
causes a blue shift of 17.4 cm~^, an order of magni- 
tude larger than what is commonly observed. Within 
the limits imposed by the coarse resolution of the spec- 
trum (~ 1 cm~^), the excited-state inversion splitting is 
found to be decreased by 31%. At the time, only specu- 
lations could be made as for the cause: it was suggested 
that helium might increase either the tunneling barrier 
and/or the effective mass of the oscillator, or even change 
the symmetry of the inversion potential, thus suppressing 
tunneling overall. With the knowledge we have gained in 
the meantime from the study of molecular rotation, we 
can reasonably say that the first two are going to be 
important while the third is probably not. Similar ob- 
servations have been made for the i^i and 1^2 modes of 
the ammonia dimer ||59| . In this system the inversion 
tunneling is already suppressed by the partner molecule, 
and the presence of the helium has, if any, a smaller ef- 
fect. For the same reason, the large blue shift observed 
for the monomer is not present here. Interchange tunnel- 
ing is instead affected by the presence of the helium: the 
splitting is reduced from 5-10 cm~^ to 2 cm~^. 

A capital step forward is the study of the HF dimer 
by Nauta and Miller [^: the use of a low power contin- 
uously tunable laser and the added diagnostics provided 
by Stark fields allowed a truly quantitative study of tun- 
neling motions. The interchange tunneling splitting, like 
in the ammonia dimer, is observed to be decreased, by an 
amount that is accurately quantified (-41%). Nauta and 
Miller propose that this is due to an increased barrier; the 
barrier increase is explained with the observation that the 
transition state is nonpolar, hence has less energy gain 
from polarization of the helium than the polar equilib- 
rium states. While it is difhcult to rule out additional 
contributions, such as those proposed for ammonia, this 
one must definitely be present, since the simple model 



used in Ref. ^ shows quantitative agreement with the 
experiment. Nauta and Miller observe that this interpre- 
tation also explains a rather unusual observation. The A 
rotational constant is increased (from 31.96 to 33.3 cm~^) 
upon solvation in helium, a fact which is indicative of a 
more linear structure. Such a structure exhibits a larger 
dipole moment, and is therefore consistent with a polar- 
ization energy gain. 

The experiment also measured the 1^2+1^5 and vi+v^ 
bands, from which the frequencies of the soft modes 
involving the van der Waals bond (2^4, intcrmolccular 
stretch, and 1/5, intermolecular bend) are reconstructed. 
Both are found to be increased, relative to the gas phase, 
by « 4% (- 133 cm"i up from 127.57, and - 185 cm'^ 
up from 178.67, respectively), which is a huge value com- 
pared to any other mode observed so far. This obser- 
vation is in line with the general expectation that soft 
modes will be more strongly affected by the hindering 
presence of the helium. 



C. ^He and mixed "He/^He droplets 

At present, rotationally resolved infrared spectra in 
pure or mixed '^He droplets have been only obtained for 
of SFe H, |l| and OCS |l|, ||. In mixed droplets 
the slightly lower zero point energy of "^He atoms favors 
their localization around the dopant molecule [ |140| . Thus 
a small doped *He cluster within the large '^He droplet 
is formed; clusters with a small number of ^He atoms 
{N4 — 20-30) produce considerably broadened spectra 
(ascribed to inhomogeneous broadening induced by the 
spread of N4). Otherwise the resulting spectra are, to 
first approximation, the same as in pure ^He, except the 
temperature is lower (0.15 K), being dictated by evapora- 
tion of "^He. As the concentration of "^He in the expanding 
gas mixture is increased from 0.1 to 4% the temperature 
rises gradually to the limit value of pure ^He droplets, 
0.38 K [|l|. 

Because of the lower number d ensity of '^He, the ma- 
trix shift of the band origin (Sec. HIE ) is also affected. 
The band origin shift of SFg has been used to esti- 
mate the number of "^He atoms in the inner cluster [|ll; 
small shifts within the R branch suggest that the rota- 
tional constant of SFg in mixed droplets with a large 
number of ^He atoms (200-500) is « 3% larger then in 
pure ''He droplets. The effect is more evident for OCS, 
where th e ro t atio nal constant drops by up to 10%, at 
Ni = 60.|l4|, |l4| 

Studies of isotopically pure '^He droplets were limited 
to OCS which has a simpler Hamiltonian and a bet- 

. Whereas in pure ^He the 



ter resolved spectrum|15 



spectrum shows the well resolved rotational structure ex- 
pected for a linear molecule, in pure '^He droplets the 
spectrum collapses into a broad single band, very similar 
to what is expected in a classical fluid. Since the binary 
interaction potentials of OCS with ■^He and ''He atoms 
are virtually identical, the difference must be due to the 



21 



different states of the droplet: superfiuid for ^He and nor- 
mal fluid for '^He. Molecules in non-superfluid "^Hc are 
rapidly exchanging their energy and angular mome ntum 
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with the medium as compared to superfiuid ^He 
Thus the ability of the dopant molecule to rotate coher- 
ently over many periods is a microscopic manifestation 
of superfluidity. Recently, the spectrum in pure "^He has 
also been fitted to an unresolved rotational structure; 
this gives a rotational constant 1.6 times smaller than in 
4He.[|l42| 

The onset of superfluidity in small *He clusters has 
been investigated by intentionally doping isotopically 
pure ^He droplet containing one OCS molecule with a 
controlled number (up to 100) of *He atoms. |l^ The 
gradual appearance of a sharp rotational structure in- 
dicates that the onset of superfluidity is a smooth pro- 
cess, and is almost complete in cluster of just GO^He 
atoms, corresponding to approximately two solvation lay- 
ers. While it is understood that preservation of rotational 
coherence depends on lack of accessible excitation of the 
helium, in turn linked to the Bos e pe r mut ation symme- 
try of the first solvation shell(s) |118, 143|, the detailed 
mechanism of this effect is still to be understood. 

Further addition of the "^He atoms do not lead to quali- 
tative change of the spectrum. The appearance of partic- 
ularly sharp lines in the P branch (< 100 MHz) has been 
attributed to cancellation of rotational and vibrational 



contributions to inhomogeneous broadening. 1 141 1 



D. Molecule/hydrogen clusters in helium droplets 

The HENDI study of van der Waals clusters contain- 
ing various isotopomers of molecular hydrogen deserves 
special attention because of the highly quantum nature 
of the latter molecule. Just like the magnesium clusters 
discussed in Section IV A |^ , the arrangement and prop- 
erties of hydrogen molecules are studied through their ef- 
fect on the infrared spectrum of a larger molecule such 
as OCS |l|, H, ||. 

The first study_of an OCS-H2 (-^HD, -D2) van der 
Waals complex 
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has indicated that the hydrogen 
molecule lies in the ring-shaped locus of the potential 
minimum around OCS, with a rather large delocaliza- 
tion out of the average OCS-hydrogen plane. It is not 
clear, however, to what degree the use of Kraitchman's 
equations (which are used to derive the position of an iso- 
topically substituted atom from the resulting cha nge i n 



the moments of inertia, assuming a rigid structure 1 144 



is justified given the large amplitude of the hydrogen vi- 
bration. However, the resulting OCS-H2 structure was 
found to be in reasonable agreement with the predictions 
of recent ab initio calculations ]14E(| . 

The same authors have found jBSf that upon addition 
of more para-hydrogen molecules to the 0CS-(pH2)„ 
van der Waals complex in pure *He (0.38 K) and mixed 
^He/'^He (0.15 K) droplets, the ro-vibrational spectrum 
rapidly approaches that of a symmetric top, until, for 



n = 5 and 6, the Q-branch disappears and the symmetry 
of the spectrum is that of a linear molecule again, just 
as for a pure OCS molecule in a helium droplet. This is 
most probably due to the fact that the bosonic (J = 0) 
pH2 molecules form a complete ring of Cnv (n — 5-6) 
symmetry around OCS. As discussed for helium excita- 
tions in Section III, aziniuthal excitations of such a ring 
are forbidden if K is not an integer multiple of n. In the 
cold environment of the helium droplets, such high-i^ 
levels are not populated, and there can be no Q-branch. 



In the same article, a very interesting comparison 
has been made with the same complex but using or- 
f/io-deuterium (0D2) instead of pH2. One would expect 
that the 6-fold nuclear spin degeneracy of 0D2 makes the 
probability for building a ring of indistinguishable 0D2 
molecules equal to 6^~", which amounts to 0.8% for n = 5 
and 0.1% for n = 6. The resulting decreases in intensity 
of the Q-branches of OCS-(oD2)„ complexes are there- 
fore negligible. Indeed, the spectra of all OCS-(oD2)„ 
clusters up to n = 8 feature a Q-branch. 



Most surprisingly, the same authors have found pl| ] 
that the Q-branch disappears again for 0CS~(pH2)„ 
clusters with n = 14-16 at 0.15 K, but not at 0.38 K. 
For these values of n, the pH2 molecules probably form 
several rings around OCS, with different numbers rii of 
molecules per ring. One could imagine a situation where 
two rings are formed, with ni — 6 and ?i2 = 8 pH2 
molecules, with a C2V structure. Based upon the re- 
sults for the OCS-H2 monomer, a moment of inertia of 
la « 370 u A2 can be estimated for this structure. This 
would imply that the lowest accessible state with rota- 
tional angular moment along the OCS axis {K — 2) will 
have an excitation energy of 0.26 K, allowing significant 
population at 0.38 K but not at 0.15 K. The next two 
pH2 molecules could then be added at both ends of the 
OCS on the molecular axis, preserving the C2„ structure, 
and finally resulting in a filled shell of 16 pH2 molecules 
around the OCS (com pare to a filled ''He shell that con- 
sists of 17 molecules 



107]). This picture neglects any 



additional reduction of the rotational constant by the 
helium, which is expected to be small because such a 
highly symmetric structure may not induce much angu- 
lar anisotropy in the helium. The authors of Ref. |5l| 
propose the radically different explanation that at the 
lower temperature (0.15 K), large delocalization of the 
pH2 molecules leads to bosonic exchange and thus to a 
"superfiuid" state of the hydrogen rings, which effectively 
makes the symmetry group Coov and forbids any states 
with K ^ 0. At higher temperature (0.38 K), however, 
delocalization is less, and thermal population of \K\ > 
states yields a Q-branch. The field of (pH2)„ clusters will 
need further study in order to corroborate this claim of 
superfiuidity. 
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V. CONCLUSION 

HENDI spectroscopy has evolved from earlier molec- 
ular beam techniques for isolated and complcxed atoms 
and molecules and has, in a rather short time, reached 
the point where it can be applied to the study of a wide 
variety of problems. Ro- vibrational spectroscopy in this 
novel extremely cold and finite environment reveals rich 
spectra, often with much more detail than in other con- 
densed phases. The field has made significant progress in 
developing the experimental and theoretical tools needed 
to extract relevant information from these spectra. How- 
ever, a number of fundamental questions remain which 
principally relate to the nature of the interactions of the 
molecular motions with the highly quantum many-body 
dynamics of the helium. We still do not understand the 
mechanism by which molecules come into equilibrium 
with the droplets, nor do we know precisely the time 
scales for many of the relaxation processes that occur. 

Despite the fact that much work remains to be done 
on these basic issues, the field is also rapidly making the 



transition from using molecules to study the physics of 
the droplets to using droplets to create and study novel 
chemical species. Vibrational spectroscopy has produced 
the bulk of the results in this direction, due to the high 
resolution available and the nearly additive effects from 
multiple perturbers. New tools, especially continuously 
tunable IR lasers at longer wavelengths than 3.5 /xm and 
the development of sources of radicals for pickup by the 
droplets, promise to greatly expand the range of Chemi- 
cal Physics that can be done with the smallest and coldest 
chemical vials yet devised. 
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